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SECTION  I 


INTRODUCTION 


The  primary  purpose  of  this  development  contract  was  an  EPIC 
link  with  HULL  to  allow  the  Ballistics  Research  Laboratory  to  run 
multiplate  penetration  problems  with  relative  ease.  The  link  is 
patterned  after  a  similar  link  between  HULL  and  TOODY  for  two- 
dimensional  problems. 

In  addition  to  the  linking  routines,  a  great  deal  of  time  was 
spent  in  making  EPIC  simpler  to  run  in  a  production  mode.  Several 
EPIC  routines  were  combined  under  the  SAIL  update  system.  As  a 
result  EPIC  routines  are  automatically  sized  to  use  the  minimum 
core  possible.  Sizing  includes  array  sizing  as  well  as  insertion 
or  deletion  of  large  blocks  of  coding.  EPIC  time  history  station 
routines  were  modified  to  allow  the  user  to  input  derived  positions 
(instead  of  node  and  element  numbers) ,  to  save  data  based  on  a  per¬ 
centage  change  in  key  variables  since  the  last  save,  and  to  pack 
the  data  into  optimum  size  buffer  areas.  The  time  history  plotter 
was  modified  to  allow  any  amount  of  data  to  be  plotted  by  deleting 
the  requirement  that  the  data  for  a  station  fit  in  memory. 

In  addition,  a  more  realistic  metal  equation  of  state  was 
included  in  EPIC  along  with  a  failure  model  based  on  principle 
tensile  strain. 

Numerous  smaller  programming  changes  were  made  -  such  as 
format  free  input,  tape  buffer  size  reductions,  debug  points  if 
desired,  etc. 

These  programming  changes,  along  with  details  of  the  link  are 
described  in  this  report. 
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SECTION  II 


HULL/EPIC  LINK 


Lagrangian  wave  propagation  codes  such  as  EPIC  have  great 
difficulty  completing  a  plate  penetration  calculation  due  to  the 
large  distortions  which  occur  as  the  target  crushes,  strains  to 
very  large  values  and  petals  or  plugs  open.  It  is  usually  impos¬ 
sible  to  compute  such  problems  correctly  or  to  maintain  a  reason¬ 
able  time  step.  However,  this  is  the  type  of  violent  flow  for 
which  Euler ian  codes  were  developed.  They  can  handle  large  de¬ 
formations  with  ease. 

On  the  other  hand,  in  a  multiplate  target  with  long  separa¬ 
tion  distances  between  plates,  an  Eulerian  mesh  would  have  to  be 
extremely  finely  zoned  to  properly  calculate  stresses  and  strains 
in  the  slowly  deforming  projectile  as  it  passes  between  plates. 

If  the  plates  were  separated  by  any  appreciable  distance,  the  long 
flight  time  between  plates  would  make  the  Eulerian  calculation  pro¬ 
hibitively  expensive.  A  Lagrangian  code  would  excel  in  this  re¬ 
gion  since  the  deformations  set  up  by  first  plate  impact  are  rela¬ 
tively  mild.  And  the  Lagrangian  code  would  handle  this  phase  ef¬ 
ficiently  since  only  the  projectile  need  be  in  the  problem. 

The  most  efficient  multiplate  impact  calculation  would  there¬ 
fore  be  one  which  uses  the  type  of  code  most  appropriate  in  each 
time  regime.  This  is  the  basic  rationale  behind  the  development 
of  the  HULL/EPIC  link. 

HULL  will  handle  the  first  plate  penetration  and  then,  in 
effect,  hand  the  problem  over  to  EPIC  for  the  long  flight  to  the 
next  plate.  Upon  reaching  this  next  plate,  EPIC  can  hand  the 
problem  back  to  HULL  for  the  plate  penetration.  If  there  is 
penetrator  left  after  this  second  plate  penetration,  the  problem 
can  be  returned  to  EPIC,  then  returned  to  HULL  for  the  next  plate, 
etc . 


Actual  implementation  of  the  link  proceeds  as  follows: 

1.  A  HULL  problem  is  established  which  models  the  first 
target  plate  and  some  or  all  of  the  projectile.  Time  history 
stations  are  located  in  the  nose  region  of  the  projectile.  The 
HULL  problem  is  then  run  until  the  first  plate  is  penetrated  and 
plate  influence  on  the  projectile  ceases  or  reaches  insignifi¬ 
cant  levels. 


2.  An  EPIC  calculation  is  then  set  up  which  describes 
the  projectile  prior  to  first  plate  impact.  The  target  is  not 
modelled  in  the  EPIC  calculation.  The  calculation  is  begun  and 
EPIC  surface  nodes  are  driven  by  velocities  from  appropriate  time 
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history  station  data  stored  on  the  tape  made  during  the  HULL 
run.  Velocities  are  used  instead  of  stresses  since  the  EPIC 
boundary  must  be  driven  with  displacement  or  one  of  its  deriva¬ 
tives.  Stress  in  the  target  next  to  the  projectile  could  be  used 
but  would  require  a  complicated  routine  applying  that  stress  to 
the  appropriate  element  surface.  After  the  EPIC  calculation 
passes  through  the  time  regime  for  which  HULL  time  history  station 
data  exists,  the  nodes  which  were  driven  by  this  data  automati¬ 
cally  revert  to  free  surface  nodes  and  the  calculation  continues 
toward  the  next  plate.  Just  prior  to  reaching  the  point  in 
space  where  the  second  plate  would  exist,  the  EPIC  calculation  is 
terminated  with  a  normal  dump  which  contains  all  EPIC  element  and 
node  information  at  termination  time. 

3.  The  EPIC  dump  is  used  as  input  to  a  new  HULL  calcu¬ 
lation  which  will  describe  the  projectile  and  the  next  target 
plate. 


4.  The  EPIC  problem  can  be  continued  through  the 
second  plate  and  beyond  by  using  second  HULL  run  station  data  to 
drive  EPIC  nodes  beginning  at  the  time  of  the  last  EPIC  dump. 

5.  Upon  reaching  the  third  plate,  the  second  EPIC  run 
can  be  terminated  and  the  data  transferred  back  into  HULL  for  a 
third  plate  impact. 

6.  This  process  can  continue  indefinitely  as  long  as 
some  of  the  penetrator  survives. 

The  amount  of  projectile  which  has  to  be  modelled  in  the  in¬ 
itial  HULL  run  depends  upon  the  material  properties  of  the  pro¬ 
jectile  and  the  time  at  which  significant  plate  influence  on  pro¬ 
jectile  motion  ceases.  Enough  of  the  projectile  should  be  model¬ 
led  to  insure  that  velocity  data  in  the  HULL  time  history  stations 
is  valid.  It  would  be  undesirable  to  have  signals  return  from  the 
aft  end  of  the  projectile  to  its  nose  during  the  HULL  run,  unless 
such  signals  would  exist  in  the  fully  modelled  projectile.  The 
time  of  plate  influence  will  change  as  projectile  velocity  changes, 
its  obliquity  changes  or  the  thickness  of  the  plate  changes.  It 
is  not  possible  therefore  to  determine  a  single  time  of  influence 
to  apply  to  all  problems.  The  time  of  influence  can  be  as  little 
as  20  to  25  microseconds  for  a  projectile  impacting  a  3/8  inch  RHA 
plate  at  2  km/sec  or  as  long  as  100  to  120  microseconds  for  the 
same  projectile  slowed  to  1  km/sec  at  an  obliquity  angle  of  65  de¬ 
grees.  Whatever  the  specific  situation  the  length  of  projectile 
modelled  should  be  the  real  length  or  one  half  the  time  of  influ¬ 
ence  multiplied  by  the  projectiles'  elastic  wave  speed,  whichever 
is  smaller.  For  example,  if  the  normally  impacting  2  km/sec  pro¬ 
jectile  is  copper  which  has  an  elastic  wave  speed  of  5  km/sec  the 
length  required  to  be  modelled  is  approximately  6  centimeters  - 
unless,  of  course,  this  exceeds  the  real  length  of  the  projectile. 
If  the  obliquily  impacting  1  km/sec  projectile  is  staballoy  with 
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an  elastic  wave  speed  of  3.2  km/sec,  we  must  model  almost  20 
centimeters  of  length  to  insure  valid  HULL  station  data. 

It  would  be  highly  desirable  to  not  require  EPIC  to  "repeat" 
the  HULL  phase.  In  theory,  this  could  be  accomplished  by  ex¬ 
tracting  projectile  information  from  a  HULL  dump  made  at  the 
time  plate  influence  ceases.  In  practice  this  requires  a  three- 
dimensional  finite  element  generator  capable  of  providing  reason¬ 
able  zoning  in  a  completely  arbitrary  region.  It  would  also 
require  that  the  projectile  in  HULL  be  very  finely  zoned  so  that 
even  low  level  stresses  are  accurately  defined  at  transfer  time. 
The  former  requirement  would  be  very  difficult  to  achieve  without 
quite  a  bit  of  manual  and  time  consuming  intervention.  The  latter 
requirement  means  that  whatever  efficiency  were  gained  by  not  re¬ 
peating  the  EPIC  calculation  would  probably  be  lost  in  the  extra 
calculational  time  required  for  a  very  finely  zoned  HULL  calcula¬ 
tion.  The  best  approach  to  avoid  "repeating"  the  plate  impact 
and  to  insure  adequate  accuracy  would  appear  to  be  a  fully  inter¬ 
actively  linked  HULL/EPIC  system. 

The  most  accurate  HULL  initial  penetration  setup  is  one  in 
which  the  penetrator  is  oriented  along  the  Z-axis.  In  this  case, 
time  history  stations  should  be  used  which  move  in  a  Lagrangian 
fashion  in  the  Z  direction  and  are  Eulerian  (i.e.,  fixed  in  space) 
in  the  other  directions.  The  difficulty  in  using  purely  Lagran¬ 
gian  stations  is  simply  that  such  stations  do  not  exactly  follow 
the  flow  and  they  can  move  so  far  as  to  become  useless.  Stations 
Lagrangian  in  the  Z  direction  should  be  uniformly  placed  within 
the  projectile's  nose  region  through  the  projectile  and  radially 
to  distances  of  two  to  three  times  the  projectile's  radius  as 
illustrated  schematically  below. 


REGION 
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This  will  insure  that  stations  always  exist  in  the  projectile 
material  and  near  the  projectile's  nose  surface.  Of  course,  if 
significant  target  impact  is  expected  at  other  locations  along 
the  projectile,  HULL  stations  should  be  located  at  those  points 
also  as  might  occur  in  a  projectile  with  a  highly  flared  aft 
body  section.  The  HULL  station  tape  is  packed  to  reduce  the 
length  required  and  station  data  is  not  collected  for  a  sta¬ 
tion  if  it  has  not  changed  by  more  than  2  percent  from  the  last 
time  it  was  saved.  Therefore,  there  is  no  overriding  reason  to 
attempt  to  severely  limit  the  number  of  HULL  stations.  A  few 
hundred  HULL  stations  can  be  used  without  causing  any  calcula- 
tional  difficulties.  Even  more  stations  could  be  used  by  re¬ 
ducing  the  amount  of  information  carried  at  each  station.  Data 
carried  in  each  station  for  a  three-dimensional  run  is: 


Density  (with  station  number  masked  into 
the  lower  bits  of  the  word) 

X  position  (with  material  data  masked  into 
the  lower  bits  of  the  word) 

X  velocity 

Y  velocity 
Z  velocity 

X  acceleration 

Y  acceleration 
Z  acceleration 


x 

x 

X 

X 

X 

X 


XX 

YY 

ZZ 

XY 

XZ 

YZ 


TOTAL  STRESS  COMPONENTS 


e 

e 

e 

e 

e 

£ 


xx 

YY 

ZZ 

XY 

XZ 

YZ 


TOTAL  STRAIN  COMPONENTS 
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Y  Position 


Z  Position 

I  Internal  energy/unit  mass  (if  NVARST=23) 

The  Hull  variable  NVARST  -  contained  in  the  HULL  Z  block  - 
indicates  the  number  of  variables  carried  per  station.  In  the 
normal  HULL  run  this  is  set  to  22.  The  user  can  reset  it  to 
23  and  HULL  will  automatically  also  carry  internal  energy  per 
unit  mass.  It  is  a  relatively  simple  matter  to  change  HULL  to 
carry  fewer  variables  -  just  the  station  positions  and  velocities 
for  instance  -  and  to  change  EPIC  to  accept  this  reduced  amount 
of  information.  NVARST  could  be  changed  to  a  smaller  number  in 
KEEL  input  and  coding  deleted  or  changed  in  HULL  subroutine 
STATON  which  causes  variables  to  be  saved  in  the  STAD  array.  The 
EPIC  changes  are  even  simpler.  In  EPIC  subroutine  RDH,  one  need 
change  only  the  arguments  in  the  buffered  in  AHULL  array  in 
loops  in  which  this  information  is  transferred  to  the  EPIC 
variables  XHULL ,  YHULL ,  ZHULL,  UHULL ,  VHULL  and  WHULL.  These 
latter  arrays  contain  HULL  station  X,  Y,  Z,  U,  V  and  W  values. 

The  material  information  placed  in  the  lower  bits  of  the  X 
position  is  binary  information  which  describes  whether  HULL  mater¬ 
ials  are  in  the  zone  (1)  or  not  in  the  zone  (0)  at  the  time  of 
station  updating.  For  example,  if  HULL  is  running  with  three 
materials  (CU,  RHA  and  AIR)  and  the  HULL  material  numbers  are: 

CU  =  1 
RHA  =  2 
AIR  =  3 

Then  if  a  zone  contains  only  air  the  lower  bits  in  the  density 
word  will  be: 

100 

If  the  zone  contains  air  and  copper  the  bits  will  be: 

101 

If  all  three  materials  are  contained  in  the  zone,  the  bits  are: 

111 

As  seen,  the  material  data  is  stored  from  right  to  left  in  bits 
which  correspond  to  the  material  number.  When  EPIC  examines  a 
HULL  station,  it  unmasks  the  material  type  word  and  decides  if 
the  word  indicates  that  HULL  projectile  material  was  contained 
in  the  zone  at  the  time  the  station  was  updated.  The  HULL 
material  numbers  for  projectile  materials  are  carried  in  the 
MATHUL  array  which  is  input  to  EPIC  at  the  time  the  EPIC  MAIN 
program  is  run.  Because  the  HULL  material  type  word  contains  15 
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binary  bits,  the  number  of  materials  identified  in  the  word  is 
limited  to  15.  In  the  case  illustrated,  if  the  projectile  is 
copper  then  MATHUL  would  be  1.  The  material  type  word  is 

logically  added  (AND)  to  2  (MATHUL_1) #  jf  the  reSult  is  greater 
than  zero  projectile  material  exists  in  the  station  and  the 
station  can  be  used  to  drive  EPIC  surface  nodes.  (Recall  that 
in  an  "AND"  operation,  1  and  0  become  0,  0  and  1  become  0,  0 
and  0  remain  0.  Only  1  and  1  will  become  1  —  i.e.,  a  number 
greater  than  zero) .  In  a  linked  run  the  ISLAVE  node  variable 
is  used  to  contain  EPIC  material  number  for  the  associated  ele¬ 
ment.  The  EPIC  MAIN  program  not  only  asks  for  a  MATHUL  array, 
but  also  for  a  corresponding  MATEPIC  array.  The  arrays  indicate 
the  relationship  between  HULL  materials  and  EPIC  materials. 

An  EPIC  node  is  driven  only  by  HULL  stations  which  contain  some 
of  the  appropriate  material.  The  relationship  in  this  array 
can  be  changed  after  the  second  and  subsequent  HULL  runs. 

HULL  station  data  is  contained  in  two  arrays.  One  array 
contains  data  valid  at  some  HULL  time,  tj,  and  one  array  contains 
the  data  at  the  next  HULL  time,  t2 ,  where  EPIC  time  is  equal  to 
one  of  these  limits  or  between  them. 

When  EPIC  requires  HULL  velocity  for  a  node  it  searches 
through  both  arrays  for  all  stations  which  contain  the  appro¬ 
priate  material  and  which  are  within  a  DELTA  distance  from  the 
node.  DELTA  is  set  to  1  centimeter  but  can  be  changed  to  any 
value  by  changing  the  data  statement  in  MAIN  subroutine  HVEL . 

Note  that  for  a  linked  run,  EPIC  must  use  the  CGS  system  of 
units.  However,  for  an  unlinked  EPIC  voyage  under  SAIL,  any 
consistent  set  of  units  is  permitted.  The  station  nearest  to 
the  node  is  used  if  it  contains  the  correct  material  at  both 
times.  There  is  no  spatial  interpolation  employed.  Once  the 
station  has  been  found,  velocity  components  are  established  for 
the  node  by  linearly  interpolating  between  data  at  ti  and  t2  to 
find  data  valid  at  the  current  EPIC  time.  This  HULL  velocity 
data  is  then  used  directly  to  drive  the  node  in  question.  If  a 
station  cannot  be  found,  the  X-component  of  velocity  is  returned 
as  -7777777  and  velocity  is  computed  for  the  node  as  would  nor¬ 
mally  be  accomplished  using  the  forces  on  the  elements  surround¬ 
ing  the  node.  When  HULL  data  ceases,  the  variable  NHULL  is  set 
to  zero  and  node  velocities  are  calculated  solely  from  forces 
on  the  elements. 

NHULL  was  set  to  1  by  the  EPIC  generator  when  the  user  asked 
that  surface  nodes  for  certain  projectile  regions  be  identified 
as  being  driven  by  HULL  input.  Nodes  on  the  surfaces  of  rods  and 
nose  sections  can  be  so  identified.  Two  generator  input  variables 
are  used  for  this  purpose  for  a  rod.  If  NHSURF  is  set  to  1,  the 
generator  assumes  that  all  cylindrical  surface  nodes  will  be 
driven  by  HULL  -  these  are  all  front  surface  nodes  including 
those  on  the  cylindrical  surface.  If  both  NHFRNT  and  NHSURF  are 
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1  all  surface  nodes  -  except  those  on  the  aft  face  -  are 
driven  by  HULL  input.  These  regions  are  illustrated  below 
where  dots  identify  nodes  driven  by  HULL  input. 


NHSURF=1 


NHFRNT=1 


NHSURF=1 

and 

NHFRNT=1 


For  nose  sections  there  is  only  one  input  variable  -  NHSURF.  If 
it  is  1  all  nose  section  surface  nodes  are  driven  by  HULL  input. 

Nodes  to  be  driven  by  HULL  stations  are  identified  with  a 
1  in  the  ten  thousandths  digit  in  the  constraint  word  for  the 
node. 


It  is  somewhat  computationally  faster  to  identify  nodes  as 
driven  by  HULL  only  if  they  can  reasonably  be  expected  to  come 
in  contact  with  the  target  plate  and  only ,  of  course ,  if  they 
are  in  or  near  a  region  covered  by  HULL  time  history  stations. 
The  DELTA  check  prevents  nodes  from  being  driven  by  HULL  if 
they  are  too  far  removed  from  HULL  stations ,  but  some  computer 
time  will  be  lost  in  the  station  search  which  is  conducted  every 
cycle  for  every  node  identified  as  being  driven  by  HULL. 

The  projectile  nose  section  should  always  be  driven  by  HULL 
input.  Portions  of  the  rod  aft  of  the  nose  section  may  or  may 
not  be  so  driven  -  depending  on  the  particular  problem. 
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Of  course,  non-surface  nodes  should  never  be  driven  by 
HULL.  If  a  nose  section  exists  on  the  projectile  then  NHFRNT 
should  be  set  to  zero  for  any  attached  rod  sections . 

It  is  always  desirable  to  include  the  complete  projectile 
in  the  EPIC  calculation.  However,  in  some  impact  situations  - 
soft  metals  at  high  velocities  -  we  can  expect  the  nose  section 
to  become  highly  deformed.  This  deformation  can  drive  the  EPIC 
time  step  to  values  too  small  to  reasonably  continue  the  problem 
If  this  occurs,  the  nose  section  and  even  part  of  the  rod  can  be 
dropped  from  the  EPIC  grid  and  the  problem  continued  using  HULL 
stations  further  back  in  the  rod.  Or,  the  problem  could  begin 
with  the  nose  section  and  some  rod  missing.  In  this  case, 

NHFRNT  should  be  set  to  1  so  that  rod  front  nodes  will  be  driven 
by  HULL.  Since  the  region  which  now  constitutes  the  front  of 
the  EPIC  projectile  is  interior  to  the  HULL  projectile  the  EPIC 
calculation  will  be  an  exact  HULL  duplicate  until  the  time  that 
the  HULL  stations  are  dropped.  At  this  time  the  EPIC  calcula¬ 
tion  will  be  less  exact  since  the  EPIC  projectile  is  missing 
some  mass.  However,  the  accuracy  of  the  EPIC  run  may  be  suf¬ 
ficient  since  much  of  this  mass  probably  suffered  large  deforma¬ 
tions  and  separated  from  the  projectile  anyway.  How  much  nose 
section  can  be  dropped  from  the  EPIC  calculation  to  increase 
its  efficiency  will  clearly  be  a  function  of  the  particular 
penetration  event  being  examined  and  the  projectile  data  desired 

Nodes  can  be  dropped  from  the  projectile  as  the  EPIC  calcu¬ 
lation  is  running  through  use  of  the  material  property  EBS,  by 
dropping  projectile  surface  elements  when  they  drive  the  time 
step  below  a  given  level  or  by  dropping  all  elements  beyond  a 
given  element  number.  When  using  a  HULL  metal  equation  of  state 
in  EPIC,  maximum  principle  tensile  strain  in  an  element  is 
compared  to  EBS.  If  it  exceeds  EBS,  the  element  variable  ICHECK 
is  set  to  -1.  There  will  then  be  no  further  calculations  per¬ 
formed  for  this  element.  At  the  time  the  element  is  "dropped" 
from  the  calculation  all  nodes  defining  the  element  are  checked 
to  see  if  any  are  driven  by  HULL  input.  If  any  are  so  driven, 
then  the  remaining  nodes  -  which  will  become  surface  nodes  when 
the  element  is  dropped  -  are  modified  so  that  they  will  be 
driven  by  HULL  input.  However,  use  of  EBS  requires  some  exper¬ 
ience  with  the  particular  class  of  problems  being  run.  Indis¬ 
criminate  use  can  cause  too  much  of  the  projectile  to  be  drop¬ 
ped.  We  recommend  that  problems  be  run  without  the  EBS  check 
if  possible  (setting  EBS  to  999  in  input) .  If  this  is  not 
possible,  EBS  should  be  set  to  as  large  a  value  as  practical  for 
concluding  the  calculation.  There  is  no  input  provision  for 
changing  EBS  after  the  projectile  has  been  generated.  However, 
it  can  be  changed  with  a  simple  one  command  insertion  in  sub¬ 
routine  RECALL  after  the  material  property  record  has  been  read. 
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EPIC  MAIN  asks  for  values  for  the  variables  NHELE  and 
DTHULL  if  the  calculation  is  a  linked  one.  These  variables 
are  used  in  deciding  if  surface  elements  should  be  dropped  from 
the  projectile.  If  a  surface  element  has  an  element  number 
greater  than  NHELE  and  it  is  driving  the  time  step  below  DTHULL, 
the  element  is  dropped  (i.e.,  ICHECK  is  set  to  -1).  In  a  linked 
run,  the  EPIC  generator  uses  the  ISLAVE  digit  in  the  node  identi¬ 
fication  word  IXYZ ,  to  identify  a  surface  node.  An  element  with 
three  such  nodes  is  a  surface  element.  When  an  element  is 
dropped,  all  nodes  for  the  element  become  surface  nodes  so  that 
surrounding  elements  will  then  become  surface  elements.  The 
NHELE  and  DTHULL  inputs  are  used  primarily  to  drop  elements 
which  are  on  the  projectile  nose  and  are  being  highly  compressed 
and  in  effect,  being  taken  from  the  projectile  and  added  to  the 
target. 

If  it  is  not  desired  to  use  this  option,  set  NHELE  to  a 
very  large  number  or  DTHULL  to  zero. 

In  a  linked  run,  MAIN  allows  the  user  to  specify  up  to 
fifty  values  of  THDROP  and  NHEDRP,  where  all  elements  with 
numbers  NHEDRP  and  beyond  are  dropped  when  problem  time  equals 
THDROP.  This  was  implemented  primarily  as  a  means  of  dropping 
entire  planes  of  elements  where  NHEDRP  is  the  first  element  in 
the  plane.  However,  it  need  not  be  used  simply  for  entire 
planes. 

When  any  element  drop  option  is  used,  associated  element 
mass  is  deleted  from  the  element's  nodes  and  node  link  indicators 
are  moved  to  the  new  exposed  surface. 

Some  mild  rezoning,  such  as  element  plane  combining,  would 
be  helpful  in  maintaining  a  large  EPIC  time  step  in  high  velocity , 
soft  metal  projectile  impacts.  Such  plane  combining  was  beyond 
the  scope  of  this  contract  and  the  contract  schedule  precluded 
its  development.  However,  it  could  be  implemented  relatively 
simply  if  the  need  arises. 

The  transfer  back  into  HULL  is  accomplished  very  easily. 

The  EPIC  dump,  valid  at  transfer  time,  is  converted  into  a  form 
usable  by  HULL  through  the  EPIC  system  program  HULL IN .  This 
new  dump  then  becomes  part  of  the  input  package  in  KEEL,  the 
HULL  generator.  Details  of  setup  of  the  HULL  problem  are  dis¬ 
cussed  in  HULL  documentation.  The  tape  conversion  routine  HULLIN 
allows  the  user  to  specify  X,  Y  and  Z  bounds  in  EPIC  and  material 
numbers  for  transfer  to  HULL.  That  is,  not  all  of  the  problem 
need  be  transferred.  This  ability  to  specify  regions  for  trans¬ 
fer  means  that  surviving  portions  of  a  projectile  can  be  trans¬ 
lated  into  material  number  or  spatial  bounds.  In  addition, 
element  time  history  points  can  be  transferred  into  HULL  and 
converted  to  HULL  time  history  points  if  desired  and  if  they  are 
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in  the  material  number  and  spatial  windows.  Latest  element 
time  history  station  positions  are  printed  so  that  the  user  can 
input  them  into  HULL.  Data  from  the  EPIC  time  history  tape 
are  not  transferred.  Elements  with  ICHECK  =  -1  are  not  trans¬ 
ferred. 

HULLIN  creates  a  tape  (Tape  15)  with  the  EPIC  data  ready 
for  use  by  KEEL.  This  tape  is  not  automatically  cataloged. 

It  is  up  to  the  user  to  save  it  by  external  catalog  or  tape 
mounting  instructions. 

After  the  second  HULL  run,  the  EPIC  run  can  be  restarted 
if  there  is  any  survivable  projectile  material  left  for  travel 
to  the  next  plate.  New  MATHUL  and  MATEPIC  arrays  can  be  input 
if  required. 

This  second  EPIC  run  can  be  transferred  back  into  HULL  for 
the  next  plate  impact. 

This  process  can  continue  as  long  as  there  is  sufficient 
projectile  material  available  for  continued  EPIC  runs. 
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SECTION  III 

RUNNING  EPIC  UNDER  SAIL 


1  o 

The  SAIL  executive  language  '  is  used  to  create  and  update 
EPIC  system  routines.  SAIL  is  a  program  which  resides  in  execu¬ 
table  form  in  the  HULLIB  library.  Source  programs  for  SAIL  are 
also  available.  SAIL  is  designed  to  run  on  CDC,  Honeywell  and 
IBM  machines  with  the  change  of  one  parameter  in  its  input  file. 

SAIL  uses  a  source  file  containing  program  card  images  in 
a  special  format  and  generates  program  files  which  are  ready  for 
compilation.  The  programs  in  this  source  file  can  be  modified 
prior  to  execution  through  commands  similar  to  the  CDC  update 
system.  In  addition,  user  generated  "options"  are  available  for 
tailoring  a  program  in  many  ways.  Options  are  defined  and  given 
default  values  at  the  time  the  SAIL  source  file  is  generated. 

These  default  values  can  be  changed  in  a  number  of  ways  when 
generating  a  file  for  compilation.  The  options  can  be  used  to 
define  the  size  of  common  arrays,  to  keep  or  delete  coding  and 
even  to  be  used  by  SAIL  "if"  commands  to  define  values  for  other 
options . 

For  example,  assume  the  option  XSIZE  has  been  defined  at 
generate  time  i.e.,  when  the  program  card  images  are  first  put 
into  a  file  on  which  SAIL  can  operate.  Assume  also  that  it  is 
given  the  value  100.  When  SAIL  is  used  to  create  a  file  for 
compilation,  it  will  replace  the  expression: 

_XSIZE_ 
or  <XSIZE> 
or  (XSIZE) 

with  the  number  100  as  long  as  the  card  on  which  the  expression 
appears  has  a  $  in  column  one.  The  $  is  a  flag  which  tells  SAIL 
to  scan  the  card  for  options  delimited  with  _,  <  > ,  or  (  ) ,  signs. 
So  the  card: 

$  COMMON  /  1  /  X (_XSIZE_) 
would  become: 


COMMON  /  1  /  X (100) 


1.  AFWL  DYT  TN  75-3,  "THE  SAIL  UPDATE  AND  EXECUTIVE  PROGRAM", 

Lewis  P.  Gaby,  Daniel  A.  Matuska 

2.  AFWL  TR  78-80,  "SAIL,  AN  AUTOMATED  APPROACH  TO  SOFTWARE 

DEVELOPMENT  AND  MANAGEMENT",  Lewis  P.  Gaby, 

David  C.  Graham,  Clifford  E.  Rhoades,  Jr.,  Jan.  1979. 
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after  SAIL  processing.  If  the  default  value  of  100  is  not  de¬ 
sired  it  can  be  changed  in  the  SAIL  processing  input  deck  by 
using  the  commands: 

SAIL  OPTIONS  XSIZE=50  ENDOPTIONS  PROGRAM  MAIN 

This  input  deck  tells  SAIL  to  process  program  MAIN  from  the  SAIL 
file  and  to  change  the  value  of  the  option  XSIZE  to  50.  The 
option  value  can  also  be  changed  by  another  program  which  gener¬ 
ates  a  file  named  INPUT2 .  This  file  contains  the  following  cards 

PROGRAM  MAIN 
OPTION  XSIZE  50 

Thus  XSIZE  can  be  changed  by  a  preprocessing  program  which,  for 
example,  reads  the  input  deck  for  MAIN  and  decides  how  big  the  X 
array  should  be. 

Options  can  be  logical  or  numerical.  They  can  be  set  as 
discussed  before  or  by 

*DEF  cards  in  the  file  SAIL  will  update. 

*DEFN  XSIZE=50  will  replace  _XSIZE_  with  a  50. 

*DEFL  XLOG=TRUE  will  change  the  value  of  the 
option  named,  for  example,  XLOG  to  TRUE. 

The  options  provide  powerful  programming  tools  when  used 
with  *KEEPTO  and  *SKIPTO  commands  in  the  file  for  SAIL  updating. 
The  command : 

*KEEPTO  *36  XSIZE  EQ  50 

will  tell  SAIL  to  keep  the  next  36  cards  for  the  compile  file  if 
XSIZE  is  equal  to  50.  The  command: 

*KEEPTO  XEND  XSIZE  GT  50 

will  tell  SAIL  to  keep  all  of  the  cards  in  the  deck  until  it 
reaches  a  card: 

* LABEL  XEND 

if  XSIZE  has  a  value  greater  than  50.  Just  as  in  FORTRAN  com¬ 
mands,  SAIL  recognizes  the  operators  EQ,  GT,  GE,  LT,  LE,  and  NE. 
Logical  options  can  be  used  in  a  similar  fashion;  for  example: 

*KEEPTO  *36  XLOG 

will  tell  SAIL  to  keep  the  next  36  cards  in  the  compile  file  if 
XLOG  is  TRUE. 
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*SKIPTO  commands  are  used  in  a  similar  fashion  but  tell 
SAIL  to  skip  coding  -  instead  of  keeping  it.  For  example: 

*SKIPTO  *36  XSIZE  NE  50  tells  SAIL  to  skip 
past  the  next  36  cards  if  XSIZE  is  not  equal 
to  50 . 

Blocks  of  coding  can  be  written  once  and  then  placed  into  a  pro¬ 
gram  at  any  number  of  points.  If  the  block  of  coding  is  identi¬ 
fied  by  any  name  with  eight  or  less  characters,  say  BLOCKl,  then 
any  place: 


* INCLUDE  BLOCKl 

appears,  SAIL  will  replace  the  card  with  the  desired  coding.  A 
block  of  coding  is  given  a  name  by  the  card: 

*PR0C  NAME 

where  NAME  could,  for  example,  be  BLOCKl.  The  coding  must  begin 
with  the  *PROC  card  and  terminate  with  the  card: 

*ENDPROC 

These  procedures,  or  PROC's,  are  particularly  useful  in  generat¬ 
ing  common  blocks  for  insertion  into  several  subroutines.  For 
example,  the  following  coding  could  be  used  to  generate  a  com¬ 
mon  block. 

*PROC  /TIME/ 

$  COMMON/TIME/X (_XSIZE_) ,  Y (_YSIZE_) 

*ENDPROC 

The  common  block  with  option  values  placed  in  it  will  appear 
anywhere  a : 


* INCLUDE/TIME/ 

card  appears.  The  //  marks  are  not  required,  but  can  be  used 
and  are  helpful  in  identification  of  a  COMMON  PROC. 

A  SAIL  file  can  have  many  programs  in  it.  Each  program 
must  begin  with  a: 

*B  NAME 

card  where  NAME  is  the  name  of  the  program.  If  NAME  is  specified 
in  SAIL  input  then  only  program  NAME  will  be  processed  for  com¬ 
pilation.  There  is  one  exception  to  this.  The  program  named 
PROLOGUE  will  be  processed  each  time  any  other  program  is  pro¬ 
cessed.  Thus  coding  common  to  many  programs  can  be  placed  in 
PROLOGUE  so  that  it  does  not  have  to  be  repeated  in  each  program. 
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There  are  several  programs  in  the  EPIC  SYSTEM  file.  They 
are  as  follows: 

PROLOGUE 

PREPSIZE 

MAINSIZE 

PREP 

MAIN 

POSTl 

POST2 

EMON 

HULL INS Z 

HULL IN 

PROLOGUE  contains  common  blocks  and  routines  used  by  several 
programs. 

PREPSIZE  is  a  program  used  to  read  PREP  input  and  set  option 
values  in  an  INPUT2  file  which  SAIL  uses  when  processing  program 
PREP.  That  is,  PREPSIZE  determines  the  size  of  certain  PREP 
arrays  and  which  subroutines  must  be  kept  for  the  desired  problem. 
PREPSIZE  keeps  the  size  of  program  PREP  to  a  minimum.  It  also 
allows  interactive  debugging  of  decks  since  it  executes  the  same 
commands  as  PREP  with  the  exception  that  node  and  element  files 
are  not  actually  created. 

MAINSIZE  reads  the  input  deck  for  program  MAIN,  the  restart 
file  made  by  program  PREP  or  a  previous  MAIN  run  and  a  HULL  input 
tape  (if  a  linked  HULL-EPIC  run) .  From  these  data,  it  determines 
the  size  of  arrays  and  sets  options  for  deleting  and  maintaining 
coding  in  program  MAIN. 

PlSIZE  reads  the  restart  tape  and  time  history  tape  (if 
used)  and  determines  the  size  of  arrays  to  be  used  in  the  plot¬ 
ting  program  POSTl. 

P2SIZE  reads  the  time  history  tape  and  sets  up  options  to 
size  the  plotting  program  P0ST2. 

PREP  is  the  EPIC  generator  program  used  when  starting  an 
EPIC  problem.  It  sets  up  nodes  and  elements  on  a  restart  tape 
used  by  MAIN. 

MAIN  is  the  EPIC  numerical  physics  program  which  solves  the 
problem  described  on  the  restart  tape. 

POSTl  is  the  program  to  make  contour  plots,  vector  velocity 
plots,  etc.,  from  data  on  the  restart  tape. 

P0ST2  plots  the  time  history  data  collected  by  MAIN  on  the 
time  history  tape. 
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EMON  is  a  program  which  can  be  run  interactively  and  allows 
two-way  communication  with  a  running  EPIC  program.  Communication 
is  provided  via  a  permanent  file. 

HULL IN  is  a  program  to  convert  an  EPIC  restart  tape  to  a 
format  useable  by  HULL  in  overlaying  an  EPIC  problem  back  into  a 
HULL  grid. 

HULL INS Z  is  a  sizing  program  for  HULL IN . 

PREP,  MAIN,  POSTl ,  and  P0ST2  have  previously  been  fully 
(3  4) 

described  '  except  for  changes  made  by  OTI. 

The  options,  their  meaning  and  default  values  are  described  below: 


OPTION 

VALUES 

MEANING  DEFAULT 

VALUE 

—  COMPUTER  RELATED  — 

COMPUTER 

-1 

EGLIN  CYBER  176 

1 

0 

EGLIN  CDC6600 

1 

BRL  CDC7600 

—  TAPE 

RELATED  — 

All  tape  options  have  the  same  meaning 

0 

means  no  file 

1 

means  disk 

7 

means  7  track  tape  length 

9 

means  9  track  tape  length 

ITAPOT 

Restart  file  always  exists 

9 

ITAPLT 

NON- ZERO 

Time  History  tape  in  use 

0 

0 

Time  History  tape  not  in  use 

TAPELIB 

0 

No  automatic  tape  Library 

0 

1 

Automatic  Tape  Library  in  use 
(not  yet  implemented) 

CATPOT 

1 

Catalog  ITAPOT  file  when  full 

1 

0 

do  not  catalog 

CATPLT 

1 

Catalog  ITAPLT  file  when  full 

1 

0 

do  not  Catalog 

3.  AFATL  TR  78-81,  "Further  Development  of  the  EPIC-3  Computer 
Program  for  Three-Dimensional  Analysis  of  Intense  Impulsive 
Loading",  G.  Johnson,  Honeywell,  Inc.,  July  1978. 

4.  "Further  Development  of  EPIC-3  for  Anisotropy,  Siding  Surfaces, 
Plotting  and  Material  Models",  G.  Johnson,  Oct.  '79  (to  be 
published  as  a  BRL  report) 
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Cataloging  for  ITAPOT  and  ITAPLT  will  occur  when  the  file 
has  reached  the  length  defined  by  ITAPOT  and  ITAPLT  or  at  the 
end  of  the  current  run  —  whichever  occurs  first.  If  the  file 
is  cataloged  because  of  its  assigned  tape  length,  a  new  file 
will  be  started  and  the  run  will  continue.  The  files  are  cata¬ 
loged  as: 

EP IC  3TAPE  9XXXXPTY YY Y 
and 

EP IC  3TAPE  7 XXXXPTY YYY 

where  the  case  number  (now  a  floating  point  variable)  for  the  pro¬ 
blem  being  run  is  XXXX.  YYYY.  The  files  are  cataloged  in  a 
directory  determined  by  the  USER  option.  At  BRL ,  they  are  cata¬ 
loged  on  the  z  machine. 

On  a  restart  the  latest  tape  cycle  is  automatically  attached 
and  used  to  begin  the  problem.  It  is  returned  after  the  problem 
begins.  So  every  restart  will  generate  at  least  one  more  cycle 
of  either  tape.  A  given  cycle  of  the  restart  file  may  contain 
more  than  one  EPIC  restart  dump.  It  is  up  to  the  user  to  insure 
that  he  does  not  violate  any  local  prohibitions  on  the  total 
number  of  cycles  cataloged. 

CATPLOT  1  Catalog  Tape99  from  PCSTl  or  P0ST2  1 

0  do  not  catalog 

The  files  are  cataloged  as: 

P 1XXXXPTYYYY 
and 

P  2  XXXXPTY YYY 

where  the  PI  file  is  from  P0ST1  and  the  P2  file  is  from  P0ST2 . 

They  are  cataloged  under  a  directory  determined  by  USER.  At  BRL 
the  files  are  cataloged  on  the  A  machine. 

USER  User  works  to  determine  directory 

Name  for  cataloging  at  BRL  and  Eglin  AFB 

-  REMOTE  CONTROL  RELATED  - 

RCONTRL  0  No  remote  control  0 

1  Remote  Control  is  On 

The  remote  control  option  is  explained  in  detail  under  the  pro¬ 
gram  EMON  description: 

-  VECTORIZED  CODING  RELATED  - 

VECTOR  0  Vectorized  coding  not  desired  0 

(only  available  for  linked  HULL  runs) 

1  Vectorized  coding  desired 
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DEBUG  RELATED 


DEBUG 

SLIDE 

NHULL 

EPICEOS 

HULLEOS 

HE 

CONCRT 

ISOTRP 

MATLIB 

PYFAIL 


0  No  debug  print  desired  0 

1  Debug  print  desired  in  whatever 

program  is  being  run 

-  PROBLEM  RELATED  - 

0  No  sliding  surface  in  the  problem  1 

1  There  is  a  sliding  surface 

0  No  HULL  input  tape  0 

1  HULL  tape  required  (Tape  14)  for 

linked  run 

0  Normal  EPIC  metal  equation  of  state  1 

(MTYPE=1  or  2)  not  used 
1  Normal  EPIC  metal  equation  of  state 

is  being  used 

0  The  HULL  equation  of  state  (MTYPE=6)  0 

is  not  being  used 

1  The  HULL  equation  of  state  is  being 

used 

0  The  EPIC  high  explosive  equation  0 

of  state  (MTYPE=3)  is  not  being 
used 

1  It  is  being  used 

0  The  EPIC  concrete/geology  equation  0 

of  state  (MTYPE=4)  is  not  being  used 
1  It  is  being  used 

1  The  EPIC  anisotropic  equation  of  state  1 

(mtype=5)  is  not  being  used 
0  It  is  being  used 

1  Obtain  material  properties  from  the  0 

material  library,  (not  yet  implemented) 

0  Input  all  material  properties 

0  Do  not  use  a  P/Y  vs  Strain  to  1 

failure  surface  for  failure  deter¬ 
mination 

1  Use  a  P/Y  surface 

(PYFAIL  can  only  be  used  with  the  HULL 
equation  of  state.  If  HULLEOS  is  1, 
PYFAIL  is  automatically  set  to  1.  If 
other  materials  are  used  in  addition  to 
those  under  the  HULL  equation  of  state, 
they  will  use  the  standard  EPIC  failure 
model) 
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STRAIN 


FAIL 


NROD 

NNOSE 

NPLATE 

NSPHER 

LROD 

LNOSE 

LPLATE 

LSPHER 

ILPLOT 

ISPLOT 

NNPLOT 


0  Do  not  save  strains  for  elements  0 

1  Save  total  strains  for  elements 

(Total  strains  are  required  for 
the  PYFAIL  option.  If  PYFAIL=1, 

STRAIN  components  are  automatically 
saved) 

0  If  failure  occurs  (ICHECK=1)  do  not  0 

change  the  element  physics 
1  If  failure  occurs,  do  not  let  the 

element  have  strength  or  a  tensile 
pressure  capability. 

•  -(FAIL  only  has  meaning  if  the  HULL 
equation  of  state  is  being  used.  The 
other  EPIC  equations  of  state  will 
automatically  change  element  physics 
if  failure  occurs) 

0  Subroutine  NROD  will  not  be  used  1 

1  Subroutine  NROD  will  be  used 

0  Subroutine  NNOSE  will  not  be  used  1 

1  Subroutine  NNOSE  will  be  used 

0  Subroutine  NPLATE  will  not  be  used  1 

1  Subroutine  NPLATE  will  be  used 

0  Subroutine  NSPHER  will  not  be  used  1 

1  Subroutine  NSPHER  will  be  used 

0  Subroutine  LROD  will  not  be  used  1 

1  Subroutine  LROD  will  be  used 

0  Subroutine  LNOSE  will  not  be  used  1 

1  Subroutine  LNOSE  will  be  used 

0  Subroutine  LPLATE  will  not  be  used  1 

1  Subroutine  LPLATE  will  be  used 

0  Subroutine  LSPHER  will  not  be  used  1 

1  Subroutine  LSPHER  will  be  used 

0  Do  not  save  internal  loads  data  on  0 

tape  ITAPLT 

1  Save  internal  loads  data 

0  Do  not  save  system  data  on  tape  ITAPLT  0 

1  Save  system  data 

The  number  of  nodes  at  which  time  0 

history  data  is  desired  (no  longer 
limited  to  20) 
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NLP LOT 

The  number  of  elements  at  which 
time  history  data  is  desired  (no 
longer  limited  to  20) 

0 

ITIME 

0 

Time  history  data  is  not  being 
saved 

0 

1 

Time  history  data  is  being  saved 

RIGID 

0 

There  is  no  rigid  body  in  the 
problem 

0 

1 

There  is  a  rigid  body  in  the  problem 

COMBINE 

0 

Do  not  combine  projectile  planes 

0 

1 

Combine  planes  (not  yet  implemented) 

LINK 

1 

The  calculation  is  (or  was)  a  linked 

0 

run 

0 

The  calculation  is  not  a  linked  run 

— 

ARRAY  SIZING  AND  STORAGE  OPTIONS  - 

NM 

Number  of  materials  in  the  problem 

10 

NPRES 

Number  of  applied  pressure  (P,T) 
pairs 

50 

NVEL 

Number  of  applied  velocity  (V,T) 
pairs 

50 

IPOY 

Number  of  (P/Y,  strain)  pairs 

10 

LBSIZE 

Number  of  elements  in  one  band 

64 

LBAND 

Number  of  element  bands  in  temp¬ 
orary  storage  at  one  time 

1 

LASIZE 

LBSIZE  *  LBAND  =  Number  of  elements 
in  storage  area  at  one  time 

LBSP 

LBSIZE  *6+3 

LBSPP 

LBSIZE  *  18  (LBSIZE  *  24  if  STRAIN=1) 

NEBLK 

Number  of  element  blocks  in  the 
problem  (each  block  is  LBSIZE  long) 

500 

LCMEBUF 

1 

Element  buffering  arrays  will  be  in 
LCM 

0 

0 

Element  buffering  arrays  will  be  in  SCM 

LCMELEM 

1 

/ELEMNT/and/VECTOR/  arrays  and 
equivalenced  arrays  will  be  in  LCM 

0 
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NBSIZE 

NBAND 

NBANDMN 

NASIZE 

NBSP 

NBSPP 

NNBLK 

LCMNBUF 

LCMNODE 


NSTN 


0  /ELEMNT/  and  /VECTOR/  arrays  and 

equivalence  arrays  will  be  in  SCM 

Number  of  nodes  in  one  block  64 

Number  of  blocks  in  core  100 

Number  of  blocks  required  in  core  0 

at  one  time  (minimum  value  of  NBAND) 

NBSIZE  &  NBAND  =  number  of  nodes  in 
temporary  storage  area  at  one  time 

NBSIZE  +  2 

NBSIZE  *  14 

Number  of  blocks  in  the  problem 

1  Node  buffering  arrays  will  be  in  LCM  0 

0  Node  buffering  arrays  will  be  in  SCM 

0  The  /NODE/  common  blocks  will  be  in  0 

SCM  (if  NBAND  is  set  to  NNBLK  all 
nodes  will  be  in  SCM) 

1  The  /NODE/  common  blocks  will  be  in 
LCM  and  NBAND  is  set  to  the  minimum 
required  to  run  the  problem 

2  The  /NODE/  common  blocks  will  be  in 
LCM  and  NASIZE  is  automatically  set 
to  NNBLK  *  NBSIZE;  that  is,  all  nodes 
will  be  in  LCM  (NBAND=NNBLK) 

The  number  of  time  history  stations 
on  the  HULL  input  tape  (if  NHULL=1) 

-  OUTPUT  OPTIONS  - 


IPRPLOT  1  A  printer  plot  of  the  geometry  1 

will  be  produced  whenever  a  dump 
is  made  (in  PREP  and  MAIN) 

0  A  printer  plot  will  not  be  produced 

The  number  of  SAIL  options  is  rather  large.  However,  most 
are  set  automatically  by  the  sizing  program  based  on  data  in 
input  decks  and  the  values  of  other  options.  The  few  options 
actually  specified  by  the  user  normally  need  be  set  only  once — 
at  the  time  PREP  is  run.  All  option  values  are  automatically 
saved  in  records  on  both  the  restart  and  time  history  tapes. 
Therfore,  unless  an  option  value  is  being  changed  it  normally 
does  not  need  to  be  specified  again  by  the  user.  If  an  option 
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is  not  specified  by  the  user  or  set  by  a  sizing  program  it  will 
automatically  assume  its  default  value.  And  the  defaults  have 
been  set  to  reasonable  values,  so  a  user  can  simply  run  without 
setting  any  options  and  be  guaranteed  a  "good"  run  insofar  as  a 
"good"  run  is  controlled  by  option  values. 

USER,  CATPOT,  CATPLT  and  DEBUG  are  exceptions  to  the  rule 
that  options  need  not  be  reset.  The  sizing  program  must  be  able 
to  find  Tape  7  and  Tape  9  to  find  the  option  table.  Since  default 
values  of  CATPOT  and  CATPLT  are  set  to  1  -  i.e.,  files  are  assumed 
to  be  cataloged  -  the  sizing  program  must  be  given  the  correct 
USER  number  if  the  files  are  indeed  catalogued.  If  the  files  are 
not  cataloged  then  CATPOT  and  CATPLT  must  be  set  to  zero  in 
each  sizing  program  SAIL  deck,  so  that  the  sizing  program  does 
not  attempt  to  attach  non-existent  files.  DEBUG  is  reset  to 
zero  by  MAINSIZE,  P1SIZE  and  P2SIZE  prior  to  reading  new  option 
values. 

We  will  describe  how  each  option  is  set  during  a  sequence 
of  PREP,  MAIN,  P0ST1,  and  P0ST2  runs  and  we  will  indicate  which 
options  can  be  changed  by  the  user.  Prior  to  this,  however,  it 
will  be  helpful  to  understand  the  control  deck  structure  for 
running  EPIC  System  programs.  The  basic  structure  is  seen  below 
for  the  Z  machine  at  BRL. 

Job  Card 

Account  Card 


ATTACH,  HULLIB,  ID=KIMSEY 


LIBRARY,  HULLIB 


ATTACH,  CHEPIC,  ID=KIMSEY.  Brings  in  LATEST  change 

deck  for  EPIC 


COPYS,  INPUT,  TAPE 5 


Puts  input  deck  on  Tape5 


REWIND,  TAPE 5 


COPYS,  A,  TAPE 10 


COPYS,  INPUT,  A 


COPYS,  CHEPIC,  A 


REWIND,  A,  TAPE 10 


REWIND,  CHEPIC,  A 
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FILE,  OLD,  RT=S ,  BT=C 


Brings  in  EPIC  System  file 


GETPF ,  OLD,  EPICSYS , 
ID=KIMSEY,  ST=MFA. 

DYTHUL ,  I=A. 


FTN ,  I=SAIL,  B=SIZE. 


Executable  SAIL  is  DYTHUL. 
This  command  updates  EPIC 
and  makes  a  file  called  SAIL 
ready  for  compilation 

Compiles  file  SAIL  to  execute 
as  program  SIZE 


ATTACH  or  REQUEST  the  restart 
(Tape  9)  and  time  history  (Tape  7) 
tapes  if  not  running  under  automatic 
cataloging. 


SIZE. 

RETURN,  SAIL,  TAPE10,  SIZE, 
COPYS,  INPUT,  B. 

COPYS,  CHEPIC,  B. 

REWIND,  B. 

RETURN,  CHEPIC. 

FILE,  0LD,  RT=S,  BT=C. 

GETPF,  OLD,  EPICSYS, 
ID=KIMSEY,  ST=MFA. 

DYTHUL,  I=B . 

FTN,  I=SAIL,  B=MAIN . 

RETURN,  SAIL,  B. 

MAIN. 

Add  any  desired  cataloging  cards 
not  using  automatic  cataloging 

*EOR 

input  data  deck 


Executes  program  SIZE 


Brings  in  primary  program 
SAIL  update  record  and  com¬ 
bines  it  with  EPIC  change  deck 


Old  was  returned  when  SAIL 
was  executed.  It  must  be 
brought  in  again 


SAIL  is  executed  to  update 
the  EPIC  System  file 

Compiles  file  SAIL 


Executes  primary  program 


here  if 
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*EOR 

SAIL  OPTIONS  USER=10 
ENDOPTIONS 


PROGRAM  SIZE 


SAIL  UPDATE  DECK.  Word  SIZE 
should  be  PREPSIZE,  MAINSIZE, 
P1SIZE ,  P2SIZE  or  HULLINSZ  as 
appropriate 


*EOR 

SAIL  PROGRAM  MAIN 


SAIL  update  deck  for  the 
primary  program.  The  word 
MAIN  should  be  PREP,  MAIN, 
POSTl ,  POST2,  or  HULLIN  as 
appropriate . 


As  seen  in  the  deck  structure,  the  only  location  for  changing 
option  values  is  in  the  SAIL  update  record  for  the  sizing  pro¬ 
gram.  Any  option  changes  in  the  primary  program  update  record 
are  ignored.  Each  sizing  program  calls  a  subroutine  CHOPT  which 
reads  TAPE10  and  changes  option  values  through  the  INPUT2  file 
it  creates.  It  should  be  noted  that  this  is  somewhat  different 
than  the  situation  which  exists  in  running  the  HULL  program. 

In  EPIC,  all  but  derived  options  (such  as  LBSP,  NBSP ,  etc.)  are 
saved  in  an  /OPTION/  common  block  which  is  put  on  Tape9  and 
Tape  7.  The  sizing  program  reads  this  common  block  and  puts  it 
on  INPUT2  after  changing  any  values  according  to  the  sizing  pro¬ 
gram  SAIL  update  record.  When  the  primary  program  is  being  built 
by  SAIL,  any  option  values  on  INPUT2  will  override  option  values 
set  in  the  SAIL  update  record  for  the  primary  program.  Since 
essentially  all  option  values  are  specified  in  INPUT2,  were  it 
not  for  CHOPT  any  options  set  in  the  sizing  or  primary  program 
SAIL  update  records  would  be,  in  effect,  ignored.  The  HULL  pro¬ 
gram  does  not  specify  all  of  its  options  in  INPUT2  and  many  op¬ 
tions  can  be  set  in  the  primary  program  SAIL  deck.  There  are 
advantages  to  both  systems.  The  advantage  in  the  system  chosen 
for  EPIC  is  deck  simplicity.  Once  options  have  been  set  in  the 
PREP  run,  they  need  never  be  set  again  with  the  exception  of  USER 
and  possibly  CATPOT  and  CATPLT ;  whereas  in  HULL  some  options  have 
to  keep  being  reset  or  they  will  revert  to  default  values. 


We  are  now  in  a  position  to  discuss  how  options  get  set, 
automatically  and  otherwise^  during  a  full  EPIC  run  sequence. 

The  user  first  runs  PREPSIZE  and  PREP.  PREP  will  expect  to 
have  USER,  LBSIZE  and  NBSIZE  specified.  (It  actually  expects  to 
have  LBAND  also  specified  but  the  default  option  of  1  is  a  good 
value  for  all  runs) .  LBSIZE  and  NBSIZE  will  assume  values  of  64 
if  not  given  other  values  in  the  PREPSIZE  SAIL  record.  PREPSIZE 
reads  the  PREP  input  deck  and  automatically  sets: 


CONCRT ,  EPICEOS,  HE,  HULLEOS ,  ISOTRP,  IPOY ,  LNOSE , 
LPLATE ,  LROD ,  LSPHER,  NNOSE,  NPLATE ,  NROD ,  NSPHER, 

NM,  NASIZE,  NBANDMN ,  NBAND ,  NEBLK,  NNBLK ,  NPRES ,  NVEL , 
NHULL ,  LINK,  NSTN ,  PYFAIL ,  STRAIN,  SLIDE,  and  RIGID. 
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(NBANDMN  is  the  minimum  required  to  run  the  problem) 

Default  values  will  be  assumed  for: 

COMPUTER,  ITAPOT,  ITAPLT ,  CATPOT,  CATPLT ,  DEBUG, 

FAIL,  LCMEBUF ,  LCMELEM,  LCUMNBUF ,  LCMNODE ,  RCONTRL , 

VECTOR,  USER  and  IPRPLOT 

unless  otherwise  specified  by  the  user  at  this  time.  Values 
specified  by  the  user  at  this  time  will  remain  with  the  pro¬ 
blem  for  all  subsequent  runs  with  the  exception  that  LCM  avail¬ 
ability  will  be  checked  by  each  sizing  program  and  LCM  options 
will  be  reset  if  there  is  not  enough  available  LCM  storage.  At 
BRL,  we  assume  363000io  words  of  LCM  availability.  Options  set 
from  the  data  input  deck  should  not  be  changed  by  the  user,  with 
the  exceptions  of  NHULL ,  COMBINE  and  NBAND  or  NASIZE.  PREPSIZE 
sets  NBAND  and  NASIZE  (=NBAND  *  NBSIZE)  to  the  smallest  values 
possible  for  running  —  i.e.,  it  minimizes  core  storage.  If  the 
user  desires  to  run  with  all  node  storage  in  LCM,  then  he  must 
set  LCMNODE  to  2  on  the  PREP  or  subsequent  MAIN  runs . 

Program  MAINSIZE  reads  the  SAIL  options  from  the  first  dump 
on  the  restart  tape  and  examines  MAIN  input  data  to  set: 

ISPLOT,  ILPLOT ,  NNPLOT ,  NLPLOT ,  ITIME ,  NPRES  (if  IPRES=2) , 

NHULL  and  NSTN. 

It  also  checks  desired  LCM  options  against  LCM  availability  and 
resets  them  if  required. 

Program  P1SIZE  reads  the  first  dump  on  the  restart  tape  to 
optain  SAIL  options  for  the  run.  It  automatically  sets  ITIME=1 
and  resets  LCM  options  if  required. 

Program  P2ISZE  reads  the  time  history  tape  to  optain  SAIL 
options  for  the  run.  It  also  sets  ITIME=1. 
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SECTION  IV 

THE  HULL  EQUATION  OF  STATE  AND  STRESS/STRAIN  RELATIONS 


A  new  metal  equation  of  state  (MTYPE=6)  was  added  to  EPIC. 

It  is  termed  "the  HULL  equation  of  state"  simply  because  it  is 
modeled  after  that  used  in  the  HULL  code.  The  use  of  this  equa¬ 
tion  of  state  will  insure  the  consistency  required  for  linked 
HULL/EPIC  runs.  In  addition,  the  equation  of  state  and  associ¬ 
ated  stress/strain  relations  are  considered  more  physically 
appropriate  than  the  other  EPIC  metal  relations.  Since  HULL  only 
runs  in  CGS  units,  an  EPIC  linked  run  must  also  employ  CGS  units. 
In  an  unlinked  EPIC  run,  this  HULL  equation  of  state  may  be  used 
with  any  consistent  set  of  units. 


Input  for  a  HULL  Metal  is  as  follows: 


s 

r 

0 

PMIN 


initial  density 

initial  bulk  sound  speed 

shock  velocity/particle  velocity  slope 
initial  Grueisen  ratio 


-  maximum  tensile  pressure  (normally  set 
to  a  very  large  negative  number) 

v  -  Poisson's  ratio 

Ql  -  artificial  viscosity  coefficient 

Q2  -  artificial  viscosity  coefficient 

FY  -  initial  yield  strength 

FU  -  ultimate  yield  strength 

EU  -  strain  at  ultimate  yield 

EMELTO  -  energy  density  at  onset  of  melt  at 
ambient  density 

-  fusion  energy  density 


EFUS 
YF1 
EF1 
YF2 
EF2 
NPOY 
(POY,EPOY) 
EBS 


Thermal  softening  curve  parameters 


-  number  of  (P/Y,  Ep)  pairs 
-  (P/Y,  E  )  pairs 

-  maximum  tensile  strain  for  dropping  an 


element  from  the  grid  (i.e. 
ICHECK  to  -1) 


setting 


Each  of  these  input  variables  and  the  equation  of  state 
logic  will  be  explained  below: 

PRESSURE,  ENERGY  AND  ARTIFICIAL  VISCOSITY  CALCULATIONS 

The  pressure,  sound  speed,  internal  energy  and  artificial 
viscosity  are  calculated  in  subroutine  PHULL.  The  calculation 
proceeds  as  follows: 
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The  Hugoniot  pressure,  P„, 

li 


is  calculated  from: 


PR  =  Cy  +  Dy2  +  Sy3 
=  Cy 


if  V  >  0 
if  y  <  0 


where  the  excess  compression,  y,  is  defined  by: 


and  C,  D  and  S  are  calculated  from: 


C  =  po  Co2 
D  =  C  (1+2 (s-1) ) 

S  =  C  (2 (s-1) +3 ( (s-1) 2) 


The  current  Gruneisen  ratio  is  computed  from: 


with  a  maximum  value  set  at  2T  for  very  highly  distended  elements. 
This  representation  of  T  is  generally  considered  more  appropriate 
than  r  =  rQ  for  metals.  In  addition,  if  r  is  a  constant  for  all 
densities,  maxima  will  exist  in  the  total  pressure  computation 
for  relatively  small  values  of  y.  P  and  T  are  then  used  to  form: 

n 


FI  =  PR  (1  -  T  y/2) 


and 


F2  =  Tp  =  Topo 


which  eventually  will  be  combined  to  form  the  Mie-Gruneisen 
total  pressure: 


P  =  FI  +  F2  •  E 


when  E  is  the  element's  internal  energy  density. 

Prior  to  combining  Fl  and  F2  to  form  the  element's  pressure, 
these  quantities  are  subjected  to  two  checks  to  keep  them  at 
reasonable  levels  for  very  highly  expanded  elements. 

If  -0.25  <  y  <  -0.2,  then  Fl  is  varied  from  its  value  at 
y  =  -0.2  to  zero  at  y  =  -0.25.  The  theory  behind  this  check  is 
simply  that  any  metal  expanded  to  a  density  of  0.8  times  initial 
density  must  have  failed  and  an  element  with  such  conditions  must 
consist  of  highly  fractured  metal  particles  with  no  overall  ten¬ 
sile  capability. 
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If  -0.4  <  y  <  -0.25,  then  F2  is  varied  from  its  calculated 
value  at  y  =  -0.25  to  zero  at  y  =  -0.4.  The  purpose  of  this  check 
is  to  prevent  extremely  expanded  elements  from  existing  at  high 
pressure  levels  due  to  large  internal  energy  values. 

These  two  checks  taken  together  keep  an  element  from  main¬ 
taining  high  tensile  or  compressive  pressures  when  it  is  in  a 
very  higly  expanded  state.  The  Fl  check  provides  the  type  of 
reasonable  tensile  pressure  cutoff  normally  provided  by  a  PMIN 
check.  Thus  PMIN  is  redundant  and  is  only  included  in  the  equa¬ 
tion  of  state  to  provide  flexibility  for  the  user. 

After  the  Fl  and  F2  checks  are  complete,  pressure  is  com¬ 
puted  from: 


P  =  Fl  +  F2  •  E 


and  a  sound  speed  squared  is  computed  from: 


SS2  =  +  P*(F2)/p2  +  (4/3)  G/ po 

dp 


where  G  is  the  shear  modulus. 

The  equation  for  SS2  takes  advantage  of  the  thermodynamic 
identity 


where  subscript  s  refers  to  entropy.  SS2  is  then  used  to  compute 
the  artificial  viscosity  from  the  standard  EPIC  relationship. 

After  these  computations,  P  is  checked  against  PMIN.  If  P 
is  less  than  PMIN,  both  P  and  E  are  reset.  If  the  SAIL  option 
FAIL  is  set  to  1,  the  equation  of  state  checks  P  and  limits  it 
to  compressive  values  if  ICHECK  =  1.  As  a  last  check  on  P,  the 
equation  of  state  compares  the  computed  P  to  a  variable  PMAX,  now 
set  at  100  megabars.  Values  of  pressure  this  large  should  not  be 
exceeded  if  the  code  is  running  stably. 


STRESS/STRAIN  CALCULATIONS 


Subroutines  involving  stress  and  strain  calculations  were 
charged  to  make  input  simpler,  to  incorporate  a  physically  realis¬ 
tic  failure  model  and  a  thermal  softening  model  consistent  with 
that  used  in  the  HULL  code. 

The  shear  modulus,  G,  is  now  computed  from: 
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where  v  is  the  input  value  of  Poisson's  ratio  and  is  the 
initial  bulk  modulus,  poCo2. 


The  non-thermally  softened  yield  strength  is  calculated  from 


y  -  py  +  (f„-fx)  • 
y  =  F0  i£  ep  >  Eu 


ep/e0  i£  ep  1  eu 


where  ep  is  the  accumulated  plastic  strain  for  the  element. 

This  value  of  Y  is  then  multiplied  by  a  thermal  softening 
factor  which  varies  from  1  when  the  internal  energy  density  is  0 
to  0  when  the  internal  energy  density  is: 


E,.  +  0.5 
M 


EFUS 


Where  E^  is  the  onset  of  melt  energy  at  current  density  given  by: 


EM  emo 


(l  +  C6  •  n  +  C7  •  n2) 


(5) 


where  E^q  is  the  energy  density  at  the  onset  of  melt  at  ambient 
density. 

C6  and  C7  are  material  coefficients  calculated  from: 

C6  =  2  To  -  2/3 

C7  =  (To  -  1/3)  (2  •  To  +  1/3)  -  1 

ri  =  volumetric  strain  =  1  -  po/p 

and  Epus  is  the  energy  density  required  for  completion  of  melt 
at  ambient  density. 

The  thermal  factor  is  described  by  means  of  a  linear  fit 
between  four  points  as  shown  below: 

YF2 

YF2- 

THERMAL  FACTOR 

a 

0 
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The  points  (EF1,  YF1)  and  EF2,  YF2)  are  input  values.  If  only 
one  point  is  desired,  both  points  should  be  input  as  the  same 
values . 

There  is  no  explicit  strain  rate  factor  used  in  the  yield 
strength  calculation.  It  is  assumed  that  the  material  will  be 
deformed  at  strain  rates  on  the  order  of  100  sec"-*-  or  greater. 
The  yield  strength  of  most  metals  does  not  change  appreciably 
with  strain  rate  beyond  this  point. 

The  failure  model  used  with  a  HULL  metal  in  EPIC  is  a 
generalized  model  which  should  correctly  predict  failure  for 
any  state  of  stress.  The  user  defines  a  table  of  P/Y  vs  plastic 
strain  at  failure  points.  The  type  of  failure  surface  generated 
should  generally  appear  as  seen  below  (if  all  states  of  stress 
are  to  be  considered) : 


P/Y 


5.  Richard  Grover,  "Liquid  Metal  Equation  of  State  Based  On 
Scaling",  Journal  of  Chemical  Physics,  Vol.  55,  No.  7, 

1  October  1971 


30 


The  curve  illustrated  above  is  developed  from  four  input  pairs. 
The  basic  theory  behind  such  a  curve  is  that  failure  will  occur 
at  small  values  of  strain  for  essentially  plane  strain  states 
(tensile  P/Y  values  of  —1  to  —2) ,  at  medium  values  of  strain  for 
plane  stress  states  (tensile  P/Y  values  from  -0.2  to  -0.4)  and 
at  very  large  values  of  strain  for  compressive  stress  states. 

The  strain  used  for  the  model  is  the  maximum  principle  plastic 
strain  if  this  strain  is  positive  -  i.e.,  if  it  is  tensile. 

The  model  requires  calculation  of  the  maximum  principle 
plastic  strain,  which  means  that  current  strains  must  be  known 
for  each  element.  If  the  HULL  equation  of  state  is  used,  the 
SAIL  option  STRAIN  is  automatically  set  to  1  and  the  isotropic 
strains  exx,  eYY'  ezz'  eXY'  £XZ'  and  eYZ  are  saved  for  each  ele¬ 
ment.  The  maximum  principle  strain  is  computed  by  finding  the 
maximum  root,  E,  of  the  determinant  equation: 


det 


this  determinant  equation  reduces  to: 


e3  +  A  •  e2  +  B  •  e+C=0 


where : 


+  P  p  2  4-  p  p  2 

YY  XZ  EZZ  eXY 

4-  C  C~  4-  C~  r*  ^ 


"eXX  ~  eYY  ”  ezz 
EXX  £ZZ+  £XX  £YY  +  £YY  £ZZ 
“£XY  "  £XZ  “  £YZ 
"£XX  £YY  £ZZ  +  £XX  £YZ 

P  P  2  4-  c  c  ^ 


+  £YY  £XZ  +  £ZZ  £XY 

—  2  p  p  p 

XZ  bYZ  £XY 


Substitution  of  e '  —  e  +  A  /  3  reduces  the  equation  to 


(e')+ae'  +b=0 


where : 


a  =  ( 3B  -  A2)  /3 


and 


b  =  (2A3  -  9AB  +  27C)  /27 
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If  a  =  0  and  b  =  0,  then  all  three  roots  are  equal  and  are 
e'  =  0 

If  a  =  0  but  b  /  0,  then  all  three  roots  of  the  equation  are  the 
same  and  are  given  by: 

e •  =  b1/3 

If  a  ^  0  the  three  roots  are: 

e\  =  F  cos  (0/3) 

e'2  =  F  cos  (0/3  +  2U/3) 

e 3  =  F  cos  (0/3  +  411/3) 

where  F=2  a/3  and  cos  0  =  -(b/2)  /  V-a  3/27 .  Of  these  three 

roots  the  first  will  always  be  the  largest  or  equal  to  the 
largest. 

Once  the  largest  root  of  e'  is  found,  the  maximum  principle 
strain  is  found  by  inverting  the  previous  substitution: 

£  =  £ '  -  A/3 

The  maximum  principle  plastic  (deviatoric)  strain  is  the 
principle  strain  reduced  by  -A/3.  Note  that  this  quantity  is  £' 
before  it  is  converted.  We  employ  maximum  plastic  strain  simply 
because  it  is  the  quantity  measured  in  tests.  Spherical  elastic 
strain  can  contribute  to  failure,  but  it  is  not  present  in  the 
sample  when  the  tests  are  complete  and  the  measurements  taken. 

If  IPNODE  is  set  to  1  in  PREP  then  both  P/Y  and  maximum  ten¬ 
sile  strain  used  for  the  failure  check  are  averages  of  the  values 
at  the  element's  nodes  and  these  nodal  values  are  themselves  aver¬ 
ages  of  the  values  in  all  of  the  elements  surrounding  each  node. 
Use  of  this  option  smooths  EPIC  data  with  the  possible  reduction 
of  peak  values. 

It  should  be  noted  that  EPIC  does  not  save  internal  energy 
density  (energy  divided  by  mass)  for  an  element.  It  saves  energy 
divided  by  original  element  volume.  Thus  in  implementing  the 
equations  presented  in  this  section,  division  of  stored  energy 
by  po  was  required  in  several  locations  to  convert  it  to  energy 
density. 
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SECTION  V 

EPIC  GENERATOR  CODING  CHANGES  AND  INPUT  INSTRUCTIONS 


In  this  section  we  will  discuss  specific  coding  changes  made 
to  PREP  to  implement  the  link  and  other  improvements.  We  will 
also  discuss  the  new  PREP  sizing  programs.  Since  the  coding 
changes  do  not  affect  the  structure  of  PREP  new  flow  charts  will 
not  be  presented. 

The  PREPSIZE  program  is  run  first.  PREPSIZE  is  really  almost 
a  version  of  PREP  but  without  any  large  arrays.  Node  and  element 
geometries  are  generated  to  the  extent  required  to  determine  para¬ 
meters  such  as  the  number  of  nodes  and  elements  in  the  problem. 

The  primary  purpose  of  this  PREP  simulation  is  the  collection  of 
data  needed  to  size  PREP  arrays,  establishment  of  values  of  NBAND , 
NNBLK,  etc.  Secondarily,  the  exercise  provides  significant  debug 
of  a  PREP  input  deck  since  it  goes  through  almost  all  of  the  PREP 
instructions.  Since  it  does  not  use  any  large  arrays  it  can  be 
run  interactively  if  desired  as  a  PREP  deck  debug  program.  Nor¬ 
mally,  of  course,  it  is  run  in  a  command  chain  which  includes  the 
PREP  program. 

Input  to  PREPSIZE  (and  all  of  the  programs  in  the  EPIC  system) 
is  now  format  free.  Thus  input  decks  can  be  created  with  a  mini¬ 
mum  of  effort  on  a  remote  terminal  and  batched  to  the  main  com¬ 
puter.  The  routine  used  is  subroutine  READIN,  contained  in 
PROLOGUE.  READIN  uses  the  HULL  library  routine  VALUE  to  crack 
the  input  file  into  separate  input  words  or  values.  Each  input 
word  is  read  with  a  separate  call  to  READIN.  READIN  always  as¬ 
sumes  the  input  variable  is  floating  point  (even  though  there  may 
not  be  a  decimal  point  in  it)  so  a  separate  statement  is  required 
to  convert  it  to  integer.  Words  in  an  input  deck  are  separated 
by  spaces  or  commas.  Input  can  be  in  any  of  the  forms  below: 

X 

X.Y 
X.EY 
X.E  +  Y 
X.E  -  Y 

If  there  is  no  decimal  -  as  in  the  first  example  -  READIN  assumes 
one  at  the  end  of  the  word.  The  number  of  characters  represented 
by  X  or  Y  is  limited  only  by  the  computer's  capabilities.  READIN 
cannot  be  used  for  alphanumeric  data.  Separate  READ  statements 
are  employed  for  such  input. 

PREPSIZE  first  prints  all  input  deck  cards  with  a  call  to  the 
PROLOGUE  routine  INOUT.  Then  the  input  deck  is  rewound  and  pro¬ 
cessing  begins.  It  reads  and  processes  all  of  the  cards  to  the 
point  that  it  can  establish  values  for  the  SAIL  options: 
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NEBLK 

NBAND 

NBANDMN 

NNBLK 

SLIDE 

NM 

IPOY 

HULLEOS 

EPICEOS 

HE 

CONCRT 

ISOTRP 

PYFAIL 

NHULL 

LINK 

NVEL 

RIGID 

The  sail  options  NBSIZE,  LBSIZE  and  LBAND  are  required  to  estab¬ 
lish  NBAND,  NNBLK,  etc.  Default  values  are  used  if  new  values 
are  not  specified  in  the  PREPSIZE  SAIL  deck.  From  the  input  deck, 
PREPSIZE  also  determines  values  for  SAIL  options: 

LROD 

LNOSE 

LSPHER 

LPLATE 

NROD 

NNOSE 

NSPHER 

NPLATE 

which  are  used  as  flags  in  PREP  construction  to  determine  if  the 
element  and  node  generation  routines  of  the  same  names  are  re¬ 
quired.  Other  options  -  such  as  USER,  LCMNODE,  etc.  are  assumed 
to  retain  default  values  unless  they  are  given  new  values  in  the 
PREPSIZE  SAIL  deck. 

All  of  the  option  values  established  by  PREPSIZE  are  written 
into  the  INPUT2  SAIL  file  so  that  they  are  available  when  SAIL 
creates  the  PREP  source  file. 

If  DEBUG  is  set  to  1  in  the  PREPSIZE  SAIL  deck,  PREPSIZE  will 
print  all  of  the  data  normally  printed  by  PREP  as  it  processes 
input.  In  addition,  it  prints  node  geometry  data  for  any  nodes 
with  non-zero  values  of  IXYZ.  IXYZ  is  a  node  word  which  contains 
spatial  restraint,  sliding  surface  and  a  node  velocity  driving 
flag.  An  extra  digit  (in  the  10  thousands  place)  was  added  to 
IXYZ  as  part  of  the  link  implementation.  This  digit,  with  the 
variable  name  IVEL,  is  1  if  the  node  is  driven  by  HULL  input,  2 
if  the  node  is  driven  by  input  velocity  cards  and  0  if  the  node 
is  not  externally  driven. 
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At  the  time  PREPSIZE  is  run,  the  user  should  at  a  minimum, 
give  some  thought  to  values  desired  for  the  SAIL  options: 

USER 

ITAPOT 

CATPOT 

ITAPLT 

CATPLT 

FAIL 

STRAIN 

LCMELEM 

LCMEBUF 

LCMNODE 

LCMNBUF 

These  options  were  discussed  in  detail  in  a  previous  section  of 
this  report.  The  default  values  of  these  variables  will  provide 
an  adequate  run  but  they  may  not  provide  the  most  efficient  or 
physically  accurate  run.  The  LCM  options  can  be  changed  at  any 
time  if  desired.  The  other  options  listed  above  should  not  be 
changed  once  they  are  set  in  the  PREPSIZE  run. 

By  reading  and  processing  the  PREP  input  deck,  PREPSIZE  is 
able  to  determine  minimum  values  for  NBAND,  NASIZE  and  other 
variables  previously  estimated  by  the  user  or  left  at  very  large 
values  which  required  PREP  to  use  excess  memory.  Through  SAIL 
and  the  PREPSIZE  and  PREP  command  chain,  these  variables  are  now 
automatically  determined  and  implemented.  PREP  then  will  run  at 
the  smallest  size  possible  for  the  problem  being  generated. 

The  variables  required  in  PREP  input  are  presented  below. 
Detailed  variable  discussions  are  not  presented  for  variables 
documented  in  previous  EPIC  reports. 
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PREP  INPUT 


DESCRIPTION  CARD 

PROBLEM  DESCRIPTION, 

8A10 

MISCELLANEOUS  CARD 

CASE,  ISAVE,  IPRINT, 
NOR,  NIR,  NPL 

ISPLIT,  IPNODE*, NIR,  NNR  IZR,  NIL, 

MATERIAL  DATA 

NUMBER  AND  TYPE 


DENSITY,  SPH  HEAT,  kl,  k2 ,  k3 ,  T ,  CL,  CQ 

SHEAR  MOD,  YIELD,  ULT,  e-ULT,  PMIN ,  CRl ,  CTl,  CT2  | 

e  -SHEAR,  e-SHEAR,  e-TOTAL,  TEMP 
v 

IF  MTYPE  =  2,  EPIC  LIQUID  METAL 


DENSITY,  SPH  HEAT,  kl,  k2 ,  k3 ,  T ,  C ,  CQ 

eT-TOTAL ,  PMIN ,  TEMP 

*NOTE:  If  IPNODE  is  set  to  1  for  a  run  involving  MTYPE=6 

material,  nodal  P/Y  and  strain  to  failure  values 
are  used  in  computing  failure. 
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IF  MTYPE  =  3,  HIGH  EXPLOSIVE 


DENSITY,  ENERGY,  DET  VEL,  CL,  Cq,  e-TOTAL 

1 

Cl,  C2 ,  C3 ,  C4 ,  C5 

IF  MTYPE  =  4,  CONCRETE/GEOLOGICAL. 

J 

DENSITY,  SPH  HEAT,  TEMP,  e  -SHEAR,  e-SHEAR,  e-TOTAL, 
n  n  v 

CL'  CQ 

1 

. SHEAR  MOD,  YIELD,  ULT ,  PMIN.  CRl .  CTl,  CT2 

Cl,  C2 ,  C3 ,  C4 ,  C5 ,  C6 ,  C7,  C8 

| C9 ,  CIO,  Cll,  C12 ,  C13 ,  C14,  C15,  C16 

IF  MTYPE  =  5,  ANISOTROPIC  MATERIAL 

DENSITY,  SPH  HEAT,  kl,  k2 ,  k3,  T,  CT ,  CQ 

1 

I 

| 

1 

1 

| YIELD,  ULT,  e-ULT,  PMIN,  TEMP,  e^-SHEAR,  e-SHEAR,  e-TOTAL  | 

|0i,  02,  03,  Ex,  Ey,  Ez 

1 G  ,  G  ,  G  ,  v  ,V  ,V 

xy  yz'  zx'  xy'  yz'  zx 

1 X-YIELD ,  X-ULT,  e-ULT,  Y-YIELD ,  Y-ULT ,  e  -ULT 

1 

| Z-YIELD ,  Z-ULT,  ez~ULT,  XY-YIELD ,  XY-ULT,  Yxy-ULT 

YZ-YIELD ,  YZ-ULT,  y  -ULT,  ZX-YIELD,  ZX-ULT,  y  -ULT 

yz  1  zx 
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IF  MTYPE  =  6,  HULL  METAL 


DENSITY,  SOUND  SPEED,  S,  T ,  PMIN,  V,  CL,  CQ 

YIELD,  ULT,  e-ULT,  E  ,  Ef  ,  YFl,  EF1 ,  YF2 ,  EF2 

NPOY ,  P/Y ( 1) ,  ef(l),  ...,  P/Y (NPOY) ,  ef(NPOY) 

£MAX 

S  =  shock  velocity/particle  velocity  slope 

v  =  Poisson's  Ratio 

Em  =  melt  energy  density 

E^us  =  fusion  energy  density 

YFl ,  EF1,  YF2 ,  EF2  =  softening  coefficients 

e.,,,,  =  max  tensile  strain  allowed  before 
MAX 

dropping  element  from  grid 

TERMINATION  CARD 
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PROJECTILE  SCALE/SHIFT/ROTATE  CARD _ 

XSCALE,  Y SCALE ,  ZSCALE,  XSHIFT,  ROTATE,  SLANT 


Node  data  cards  for  projectile  (description  follows) 


TERMINATION  CARD 
0 


TARGET  SCALE/SHIFT/ROTATE  CARD 


XSCALE,  YSCALE ,  ZSCALE,  XSHIFT,  ZSHIFT,  ROTATE,  SLANT 


Node  data  cards  for  Target  (description  follows) 


TERMINATION  CARD 
0 


Element  data  cards  for  projectile  (description  follows) 


TERMINATION  CARD 
0 


Element  data  cards  for  Target  (description  follows) 


TERMINATION  CARD 

0~ 
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SLIDING  SURFACE  CARDS: 


MGEOM ,  ISEEK,  IT,  IM1 ,  IS1,  ISN,  NSLV,  NSG ,  ISR,  FRICT, 
VREF,  DREF 

MASTER  CARD  FOR  PLATE  GEOMETRY  (MGEOM=l) 

NML,  NMW,  IDL,  IDW,  IDIA 

MASTER  CARD  FOR  NOSE-ROD  GEOMETRY  (MGEOM=2) 

NOR,  NIR,  NPL 

MASTER  CARD  FOR  DISK  GEOMETRY  (MGEOM=3) 

NRING,  MCODE 

MASTER  CARD  FOR  CYLINDER  GEOMETRY  (MGEOM=4 ) 

|  NRING,  NPL,  IDL,  MCODE 

MASTER  CARD  FOR  ERODING  DISK  GEOMETRY  (MGEOM=5 ) 

NRING,  NSURF ,  IDT,  ELEl ,  "e-max,  e-ave,  AT  CHECK 

INDIVIDUAL  SLAVE  NODE  CARD(S)  -  FOR  NSG=0 

IS1,  IS2 ,  ...,  ISN 

NSG  GROUPED  SLAVE  NODE  CARD  (S) 

IS1G ,  ISNG,  INC 

TERMINATION  CARD 

0 

DETONATION  CARD 


XDET ,  YDET ,  ZDET 


INITIAL  VELOCITY  CARD 


PXDOT ,  PYDOT ,  PZDOT ,  TXDOT ,  TYDOT ,  TZDOT ,  DTl 
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NODE  INPUT  DATA 


LINE  OF  NODES  CARD 


1,  NNODE , 

XI,  Yl,  Zl,  XN,  YN,  ZN,  EXPAND, 

IX,  IY, 

IZ 

ROD  NODE 

CARDS 

IDENTIFICATION 

2 

ROD  NODE 

DESCRIPTION  CARD 

NOR,  NIR, 

NPLN ,  IRAD,  ZTOP ,  ZBOT,  EXPAND, 

NHSURF , 

NHFRNT 

ROD  NODE 

TOP  RADII  CARD  FOR  IRAD=0 

ROTOP,  RITOP,  ROBOT,  RIBOT 

ROD  NODE 

TOP  RADDI  CARD ( S )  FOR  IRAD=1 

RT  (NIR)  , 

. . . ,  RT (NOR) 

ROD  NODE 

BOTTOM  RADII  CARD(S)  FOR  IRAD=1 

RB  (NIR)  , 

.  .  .  ,  RB  (NOR) 

ROD  NODE 

TOP  SURFACE  CARD  FOR  IRAD=1 

ZT  (NIR)  , 

...,  ZT (NOR) 

ROD  NODE 

BOTTOM  SURFACE  CARD  FOR  IRAD=1 

ZB (NIR) , 

...,  ZB (NOR) 
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NHSURF 


NHFRNT 


=  0  means  nodes  on  the  surface  of  the  cylinder 
will  not  be  driven  by  input  velocity 

=  1  means  that  all  cylindrical  surface  nodes 
will  be  driven  by  HULL  input 

=  2  means  that  all  cylindrical  surface  nodes 
will  be  driven  by  velocity  input 


=  0  means  that  no  front  surface  nodes  will  be 
driven  by  input  velocity 

=  1  means  that  all  front  surface  nodes  (including 
those  on  the  front  face  at  the  outer  cylindrical 
surface)  will  be  driven  by  HULL  input 

=  2  means  that  all  front  surface  nodes  will  be 
driven  by  an  input  velocity 
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NOSE  NODE  CARDS 


IDENTIFICATION 


3 

DESCRIPTION 

NOR,  NIR, 

INOSE,  IRAD,  ROTOP,  RITOP ,  ZTOP .  ZMIN.  NHSURF 

NOSE  NODE 

TOP  RADII  CARD ( S )  FOR  IRAD=1 

RT  (NIR)  ,  , 

. RT (NOR) 

NOSE  NODE 

ZMIN  CARDS  FOR  IRftp=; 

ZM (NIR) ,  . 

,  ZM (NOR) 

NOTE:  if  NIR=0,  begin  top  radii  and  ZMIN  cards  with  RT(1) 

and  ZM ( 1 ) 


NHSURF  -  0  means  that  no  surface  nodes  are  driven 

by  external  velocities 

=  1  means  that  all  surface  nodes  are  driven 
by  HULL  input 

=  2  means  that  all  surface  nodes  are  driven 
by  a  velocity  input 
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FLAT  PLATE  NODE  CARDS 


IDENTIFICATION 


4,  TYPE 


DESCRIPTION  IF  TYPE=1 


NX,  NY,  NZ,  NXEND ,  NYEND,  IY,  X-EXPAND,  X-PART,  Y-EXPAND, 
Y-PART,  Z  ■‘•EXPAND  


DESCRIPTION  IF  TYPF.=  7. _ 

NX,  NY,  NZ,  NXEND,  NZEND ,  IY,  X-EXPAND,  X-PART,  Z-EXPAND, 
Z-PART,  Y-EXPAND 

NODE  SIZE  CARD 

XI,  Yl,  Zl,  XN,  YN  ZN 


SPHERE  NODE  CARDS 

IDENTIFICATION 

| 

1  5 _ _ _ J 

DESCRIPTION 

| 

| NOR,  NIR,  IRAD,  RO ,  RI ,  ZCG 

I 

NODE  RADII  FOR  IRAD=1 

I 

j  R  (NIR)  ,  .  .  .  ,  R  (NOR) 

NOTE:  if  NIR=0,  begin  card  with  R(l) 
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ELEMENT  INPUT  DATA 


COMPOSITE  ELEMENT  CARD 

1,  NCOMP ,  MATL,  Nl,  N2,  N3,  N4,.N5,  N6 ,  N7 ,  N8,  INC 


ROD  ELEMENT  CARDS 


IDENTIFICATION 


h 

!□ 

DESCRIPTION 

|  NOER,  NIER,  NLAY ,  Nl,  IMAT , 

MATL 

MATERIAL  IDENTIFICATION  FOR 

IMAT=1 

M (NIER) ,  . . . ,  M (NOER) 

□ 

NOSE  ELEMENT  CARDS 


IDENTIFICATION 

| 

La _ 1 

DESCRIPTION 

| 

|  NOER,  NIER,  Nl,  IMAT,  MATL  | 

MATERIAL  IDENTIFICATION  FOR  IMAT=1 

M  (NIER)  ,  .  .  .  ,  M  (NOER) 

FLAT  PLATE  ELEMENT  CARDS 
IDENTIFICATION _ 

|  4,  TYPE 

DESCRIPTION _ 

I  NLX ,  NLY ,  NLZ ,  Nl,  MATL 
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SPHERE  ELEMENT  CARDS 


IDENTIFICATION 

5 


DESCRIPTION 

NOER,  NIER,  Nl,  IMAT,  MATL 


MATERIAL  IDENTIFICATION  FOR  IMAT=1 
M  (NIER)  ,  .  .  .  ,  M  (NOER) 
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SECTION  VI 

EPIC  MAIN  CODING  CHANGES  AND  INPUT  INSTRUCTIONS 


Program  MAINSIZE  is  run  before  MAIN.  MAINSIZE  reads  the  in¬ 
put  deck  the  header  blocks  from  the  restart  tape  and  the  HULL  in¬ 
put  tape  (if  there  is  one)  and  sets  up  SAIL  options  for  processing 
MAIN.  It  also  checks  desired  LCM  options  against  availability 
and  resets  them  if  there  is  not  enough  LCM. 

Program  MAIN  has  been  modified  extensively.  Incorporation 
under  SAIL  has  been  previously  discussed  as  has  the  HULL  metal 
equation  of  state.  Of  course,  input  is  in  field  free  format. 

The  major  changes  to  MAIN  will  be  discussed  in  detail  below. 

A.  SPECIFYING  DUMP  TIME 

Input  specifying  the  time  at\ihich  restart  dumps  are  to  be 
made  has  been  simplified.  The  method  of  providing  a  card  for 
each  dump  can  still  be  used.  However,  if  only  a  single  card  is 
provided,  the  problem  time  specified  on  that  card  is  treated  as  a 
time  increment.  Of  course,  problem  termination  always  results  in 
a  dump  unless  ISAVE  is  set  to  zero. 

B.  PROBLEM  TERMINATION 

Problem  termination  can  now  occur  in  several  ways.  In  addi¬ 
tion  to  problem  time  exceeding  TMAX  or  CP  time  exceeding  CPMAX, 
or  DT  exceeding  DTMIN  a  run  can  be  terminated  by  turning  on  Sense 
Switch  1,  through  CP  time  about  to  exceed  job  card  time  limit, 
through  the  Monitor  routine  or  when  the  minimum  Z  coordinate  of 
the  projectile  reaches  an  input  ZSTOP  value  (for  linked  run  only) . 

C.  RESTART  FILE  FORMAT 

The  restart  file  format  has  been  changed  slightly.  New  HULL 
link,  equation  of  state  and  SAIL  variables  are  now  dumped.  Users 
requiring  details  of  each  record  on  the  tape  can  find  this  infor¬ 
mation  in  Subroutine  SAVE  in  the  PROLOGUE  program. 

D .  CASE  NUMBER 

Case  or  problem  number  is  now  a  floating  point  variable. 

This  provides  the  user  with  added  flexibility  in  problem  identi¬ 
fication.  Case  number  is  kept  on  the  restart  tape  and  any  time 
history  tapes.  A  comparison  of  requested  case  number  is  made 
with  the  case  number  on  such  files  prior  to  allowing  a  restart 
or  plotting  to  occur. 

E.  PROBLEM  DESCRIPTION 

The  problem  description  array  (DESC)  is  now  kept  on  the  re¬ 
start  and  time  history  tapes.  The  user  is  therefore  not  required 
to  re-enter  it  when  restarting  or  plotting  a  problem. 
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F.  TIME  HISTORY  DATA 

MAIN  has  been  rewritten  to  accept  (X,  Y,  Z)  positions  for 
element  and  node  time  plots,  to  choose  its  own  time  increment 
for  these  plots  and  to  pack  data  into  fewer  records  eliminating 
many  lengthy  end-of-record  areas. 

If  an  element  or  node  for  time  plotting  is  input  as  a  zero, 
the  code  expects  X,  Y,  Z  values  to  follow.  After  all  time  his¬ 
tory  input  has  been  read,  all  zero  node  and  element  inputs  are 
filled  in  by  a  search  through  the  grid.  Node  numbers  are  deter¬ 
mined  by  finding  the  node  at  the  minimum  distance  from  the  desired 
(X,  Y,  Z)  point.  Element  numbers  are  determined  by  finding  the 
element  within  which  the  desired  (X,  Y,  Z)  points  reside.  This 
is  accomplished  very  simply  by  calculating  the  volumes  of  the  four 
tetrahedra  defined  by  the  point  and  the  element's  four  planes.  If 
these  volumes  add  up  to  the  volume  of  the  element  then  the  point 
is  within  or  on  the  surface  of  the  element.  After  all  node  and 
element  numbers  have  been  found,  they  are  ordered  before  being 
placed  in  the  NPLOT  and  LPLOT  arrays.  These  arrays  are  now  sized 
by  SAIL  based  on  the  requested  input  so  there  is  no  artifical 
limit  imposed  on  how  many  time  history  points  can  be  specified. 

The  logic  for  determining  time  history  node  and  element  numbers 
is  found  in  subroutine  TOUT  in  PROLOGUE. 

The  user  no  longer  inputs  a  time  increment  at  which  node  and 
element  time  history  data  should  be  saved.  Node  data  is  saved  any 
time  that  the  node's  velocity  has  changed  by  two  percent  from  the 
last  saved  velocity.  Element  data  is  saved  whenever  the  element's 
mean  stress  has  changed  by  two  percent.  This  percentage  is  set  in 
data  statements  in  subroutines  LPLOTS  and  NPLOTS  and  can  be  ad¬ 
justed  to  provide  more  or  less  data.  This  technique  insures  a 
minimum  amount  of  data  to  define  all  details  of  a  complete  time 
history  trace  and  is  based  on  the  time  history  technique  employed 
in  HULL. 

Data  is  placed  in  records  approximately  1,000  words  in  length. 
When  the  array  ATIME  (1024)  would  be  over-filled  by  the  next  saved 
point,  it  is  dumped  to  tape  and  the  code  begins  to  fill  it  again. 
This  logic  is  also  contained  in  subroutine  TOUT. 

Data  saved  is  the  same  as  in  the  previous  version  of  EPIC3 
with  the  exception  that  temperature  is  always  zero  for  a  HULL 
metal  (MTYPE=6 )  and  strains  are  saved  for  an  element  if  the  SAIL 
option  STRAIN  is  greater  than  zero. 

If  the  SAIL  option  CATPLT  is  greater  than  zero,  the  time 
history  file  will  be  automatically  cataloged.  Any  previous  file 
for  the  problem  will  be  copied  to  a  new  file  at  restart  time  and 
the  new  file  will  be  extended  as  new  data  is  accumulated.  It  is 
up  to  the  user  to  purge  this  old  file  after  the  new  one  is  cata¬ 
loged  . 
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G.  HULL  LINK  ROUTINES 


Subroutine  MOTION  was  changed  to  accomodate  an  input  velocity 
or  HULL  input  tape.  If  NHULL  is  greater  than  zero  and  IXYZ  for 
the  node  is  such  that  the  node  should  be  driven  by  velocity  a  call 
is  made  to  subroutine .HVEL  to  find  node  velocities  UHULL,  VHULL, 
and  WHULL  (X,  Y,  and  Z) .  If  velocities  cannot  be  found,  the  node 
is  processed  as  a  normal  node  and  velocities  are  determined  from 
applied  forces.  If  velocities  are  found  these  are  placed  in  ap¬ 
propriate  node  storage  and  the  applied  force  calculation  is  avoid¬ 
ed.  If  NVEL  is  greater  than  zero,  IXYZ  is  checked  and  the  input 
velocity  is  applied  if  appropriate. 

Subroutine  HVEL  directs  all  HULL  input  tape  activities.  When 
it  is  first  called,  HVEL  calls  RDH  to  read  the  header  blocks  from 
the  HULL  input  tape,  Tape  14,  and  to  fill  in  initial  station 
positions  for  the  NSTN  HULL  stations.  In  addition,  the  arrays 
UHULL,  VHULL,  WHULL,  XHULL ,  YHULL  and  ZHULL  are  initialized. 

These  arrays  are  doubly  dimensioned.  The  first  dimension  corres¬ 
ponds  to  the  HULL  station  number.  The  second  dimension  is  1  or  2 
and  is  flip-flopped  as  new  data  is  read  into  core.  The  flip- 
flop  is  required  since,  in  general,  the  EPIC  problem  time  will 
be  between  two  HULL  times  at  which  data  was  saved.  Linear  inter¬ 
polation  will  be  employed  to  find  HULL  velocity  values  at  the  re¬ 
quired  EPIC  time. 

On  subsequent  calls,  HVEL  will  replenish  the  HULL  data  ar¬ 
rays  as  EPIC  time  exceeds  the  time  at  which  HULL  data  in  core  is 
valid. 

HULL  station  position  is  interpolated  between  the  two  HULL 
times  and  this  interpolated  position  is  compared  in  HVEL  to  the 
position  of  the  node  requiring  velocity  data.  Velocity  data  is 
acquired  from  a  HULL  station  if: 

a.  the  station  contains  the  projectile  material  at  both 
HULL  times,  and 

b.  it  is  the  nearest  station  containing  appropriate  material 
to  the  node  and  is  within  a  DELTA  distance  of  the  node 
(DELTA  is  set  in  a  DATA  statement  in  HVEL  to  1  centimeter) . 

If  a  station  cannot  be  assigned,  X  velocity  is  returned  as 
-7777777  and  the  node  velocity  is  determined  from  applied  forces. 

The  user  should  be  advised  that  if  the  SAIL  option  DEBUG  is 
greater  than  zero  a  great  quantity  of  data  will  output  from  this 
section  of  EPIC. 

H.  OTHER  MAIN  CHANGES 

The  HULL  metal  equation  of  state  is  described  in  another  sec¬ 
tion  of  this  report.  It  is  implemented  in  subroutine  PHULL. 
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Subroutine  BREAK  was  changed  to  employ  maximum  principle  strain 
(if  STRAIN  is  greater  than  zero)  and  a  P/Y  failure  model  to 
determine  if  an  element  has  failed.  In  addition,  it  was  modified 
to  provide  a  printout  of  pertinent  element  information  for  the 
first  100  failures.  The  ICHECK  element  word  carries  information 
concerning  element  failure.  If  ICHECK=0,  the  element  has  never 
failed.  If  ICHECK  is  greater  than  zero  it  has  failed.  If  ICHECK 
is  less  than  zero  the  maximum  tensile  strain  has  exceeded  the  HULL 
metal  input  EBS  and  is  not  subject  to  further  calculations. 

If  the  SAIL  option  FAIL  is  greater  than  zero  and  ICHECK  is 
greater  than  zero,  the  element  is  not  allowed  to  have  a  tensile 
pressure  or  any  yield  strength. 

Subroutines  STRAIN  and  ISTRESS  were  modified  to  calculate 
total  strains  from  the  incremental  strain  rates  and  the  EPIC  time 
increment  if  STRAIN  is  greater  than  zero. 

The  main  routine  was  modified  to  interface  with  the  EPIC 
Monitor  program. 

In  addition  to  these  changes,  MAIN  was  modified  extensively 
to  incorporate  SAIL  updating  features  to  reduce  core  size,  reduce 
the  number  of  "IF"  statements,  implement  desired  LCM  options,  etc. 
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MAIN  INPUT  SUMMARY 


I.  GENERAL  CONTROL  CARD 


CASE,  NCYCLE ,  IPRES,  CPMAX,  EMAX,  DTMAX , 
DTMIN ,  SSF ,  TMAX ,  NHULL 


II •  LINKED  RUN  INPUT  -  IF  NHULL  WAS  EVER  SET  TO  1 


NMLINK,  MATHUL ( 1 ) ,  MATEPIC(l), 
MATHUL (NMLINK) ,  MATEPIC (NMLINK) 


j  ZSTOP ,  VHIMP ,  THSTRT ,  NHELE ,  DTHULL 


NHDROP ,  THDROP(l),  NHEDROP ( 1) ,  .  .., 

THDROP (NHDROP) ,  NHEDRP (NHDROP) 


III. 


PRESSURE  INPUT  CARDS  -  IF  IPRES=2 


ELEl ,  ELEN,  ELEINC,  Nl,  NN,  NODEINC,  PRES 


0 


PTIME,  P(T) 


0 
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IV.  TIME  HISTORY  PLOT  CARDS 


ISPLOT,  ILPLOT,  NNPLOT,  NLPLOT ,  DTSYS 

ONLY  IF  NPLOT ( 1) =0 

NPLOT  ( 1 )  ,  X(l)  ,  Y(l)  ,  Z  (1) . 

ONLY  IF  NPLOT (NNPLOT) =0 

.  .  .  NPLOT (NNPLOT)  X (NNPLOT)  Y (NNPLOT)  Z (NNPLOT) 

ONLY  IF  LPLOT (1) =0 

1 

|LPLOT(l),  X(l)  ,  Y  ( 1 )  ,  Z(l) .  J 

| 

ONLY  IF  LPLOT (NLPLOT) =0 

| 

|  .  .  .  LPLOT (NLPLOT) ,  X) NLPLOT),  Y (NLPLOT) ,  Z (NLPLOT) 

1 

V.  DATA  OUTPUT  CARDS 


TWRITE,  ECHECK,  I LOAD ,  ISAVE 


GENERAL  CONTROL  CARD 

CASE  is  the  floating  point  problem  number. 

NCYCLE  is  the  desired  cycle  number  for  restart.  A  value 

of  -1  means  a  restart  from  the  latest  cycle. 

IPRES  controls  applied  pressure  input. 

IPRES=0  means  there  will  be  no  applied  pressure  cards. 

IPRES=1  uses  the  applied  pressure  from  the  previous 
run. 

IPRES=2  will  read  new  applied  pressure  data. 

CPMAX  is  the  maximum  CP  time  to  be  used  for  the  problem. 

This  is  a  redundant  input  now  since  EPIC  uses  CP 
time  on  the  job  card  to  terminate.  Setting  CPMAX 
to  zero  bypasses  this  feature. 
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EMAX 

DTMAX 

DTMIN 

SSF 

TMAX 

NHULL 

LINKED  RUN 
DTHULL 

NHELE 

THSTRT 

NMLINK 

MATHUL ( I ) 
MATEPIC (I) 

ZSTOP 

VHIMP 

NHDROP 

THDROP (I) 

NHEDRP ( I ) 


is  the  upper  limit  for  total  kinetic  energy 
if  applied  pressures  are  included  (IPRES=1  or 
2) .  This  input  is  used  as  a  check  for  numeri¬ 
cal  instability  and  the  code  is  stopped  if 
kinetic  energy  exceeds  EMAX.  Unlike  previous 
versions,  EMAX  is  always  input  whether  needed 
or  not.  If  set  to  zero  all  checks  will  be 
bypassed. 

is  the  maximum  integration  time  increment. 

is  the  minimum  time  increment  allowed.  If 
exceeded,  the  run  will  terminate. 

is  a  multiplier  for  the  time  step  calculation. 
It  should  be  less  than  1. 

is  the  maximum  problem  time  for  this  run. 

is  set  to  1  if  this  is  to  be  a  linked  HULL/ 
EPIC  run.  Set  it  to  0  otherwise. 


INPUT 


Minimum  time  step  allowed  for  surface  projec¬ 
tile  elements  prior  to  dropping  them. 

Elements  to  be  dropped  through  the  use  of 
DTHULL  must  be  greater  than  this  number. 

Time  to  begin  dropping  NHELE  elements  if 
DTHULL  criterion  is  met. 

The  number  of  EPIC  projectile  materials  which 
can  be  linked  to  HULL  input. 

These  arrays  define  the  relationships  between 
HULL  material  numbers  and  EPIC  material  num¬ 
bers. 

A  linked  run  stops  when  minimum  projectile  Z 
reaches  ZSTOP  assuming  an  added  Z  value  of 
VHIMP  *  TIME. 

Velocity  of  projectile  used  to  determine  if 
projectile  is  at  a  Z  value  less  than  ZSTOP. 

Number  of  (THDROP,  NHEDRP)  pairs. 

Elements  greater  than  or  equal  to  element  num¬ 
ber  NHEDRP (I)  are  dropped  when  TIME  reaches 
THDROP (I) . 
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PRESSURE  INPUT  CARDS 


SPATIAL  DESCRIPTION  CARD(S) 


ELEl  is  the  first  element  in  a  series  to  which  pres¬ 

sure  is  applied.  Must  not  be  less  than  ELE1 
or  ELEN  from  a  previous  pressure  card. 

ELEN  is  the  last  element  in  the  series  of  elements. 

It  cannot  be  less,  than  ELEl . 


ELEN INC 


is  the  increment  between  ELEl  and  ELEN. 


Nl 

NN 


is  the  node  number  opposite  the  triangular  face 
of  element  ELEl  to  which  pressure  is  applied. 

is  the  node  number  opposite  the  triangular 
face  of  element  ELEN. 


NODEINC  is  the  node  increment  between  Nl  and  NN. 

PRES  is  the  pressure  applied  to  the  triangular  faces 

of  the  elements  described  on  the  card. 

The  last  card  in  this  series  should  be  a  crad  with  a  single  0  on 
it. 


PRESSURE  MULTIPLIER  CARDS 

PTIME  is  the  time  corresponding  to  P(T) . 

P (T)  is  the  factor  by  which  all  pressures  are  multi¬ 

plied  at  the  corresponding  time. 

Cards  must  be  input  in  order  of  increasing  time. 

The  last  card  should  be  a  card  with  a  single  0  on  it. 


TIME  HISTORY  PLOT  CARDS 


GENERAL  CARD 


ISPLOT 

is 

be 

set  to 
saved. 

1  if  system  time  history 
It  is  zero  otherwise. 

data 

should 

ILPLOT 

is 

be 

set  to 
saved. 

1  if  internal  loads  time 
It  is  zero  otherwise. 

data 

should 

NNPLOT 

is  set  to  the  number  of  nodes  at  which 
data  is  desired. 

time 
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DTSYS 

is  the  time  increment  for  saving  system  and 
internal  loads  data  (if  DTSYS  is  zero  it  is 
set  to  TMAX/25) 

NODE  PLOT  CARD 


NPLOT(I) 

is  the  node  number  of  the  Ith  node  at  which 
time  history  data  is  desired. 

(X  (I)  ,  Y  (I)  ,  Z  (I)  ) 

is  the  spatial  point  at  which  node  data  is 
desired  (not  required  unless  NPLOT(I)  is  zero) 

ELEMENT  PLOT  CARD 


LPLOT  (I) 

is  the  element  number  of  the  Ith  element  at 
which  time  history  data  is  desired. 

(X(I)  ,  Y  ( I)  ,  Z  (I)  ) 

is  the  spatial  point  at  which  element  data  is 
desired  (not  required  unless  LPLOT (I)  is  zero) 

DATA  OUTPUT  CARDS 


TWRITE 

is  the  time  increment  at  which  output  will  be 
provided  (the  problem  time  of  the  restart 
dump  is  added  to  this  value  to  determine  the 
time  of  the  first  data  output) . 

ECHECK 

is  a  code  to  determine  the  type  of  printed 
output. 

If  ECHECK  is  less  than  900  system  data,  in¬ 
dividual  node  data  and  data  for  elements  with 
a  plastic  strain  equal  to  or  greater  than 
ECHECK  will  be  printed. 

If  ECHECK  is  999,  system  data  will  be  printed 
and  nodal  data  on  the  plane  of  symmetry  at 

Y=0.  No  element  data  is  printed. 

If  ECHECK  is  greater  than  1,000,  only  system 
data  is  printed. 

I  LOAD 

is  a  print  code  for  internal  loads  data. 

If  ILOAD=0  no  internal  loads  data  is  printed. 
If  ILOAD=l  internal  loads  are  printed  for  the 
cylindrical  projectile  based  on  values  of  NIL, 
NOR,  NIR  and  NPL  as  defined  in  PREP  input. 
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I  SAVE 


is  a  flag  to  determine  if  a  restart  file 
should  be  made  at  this  time  for  subsequent 
plotting. 

If  ISAVE=0,  no  file  will  be  made. 

If  I S AVE= 1 ,  a  file  will  be  made. 


There  can  be  any  number  of  these  data  output  cards.  Subroutine 
READIN  will  discover  an  end-of-file  and  not  attempt  to  read  addi¬ 
tional  cards.  In  this  case,  the  last  value  of  TWRITE  will  be 
used  as  a  continuing  time  increment  for  data  output. 
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SECTION  VII 

CREATING  A  DUMP  FOR  TRANSFER  TO  HULL 


Program  HULLIN  and  its  sizing  program  HULLINSZ  are  used  to 
create  a  Tape  15  which  KEEL  can  use  to  transfer  the  EPIC  projec¬ 
tile  back  into  HULL.  The  user  specifies  a  desired  cycle  for 
transfer  along  with  desired  materials  and  desired  spatial  regions. 
Program  HULLIN  reads  the  EPIC  dump  at  the  required  cycle  and 
creates  a  tape  for  use  by  KEEL.  Data  on  the  tape  is  as  follows: 

RECORD  DATA 

1  CASE,  NHELE ,  NHELER,  STRAIN,  XHMIN ,  XHMAX , 

YHMIN ,  YHMAX,  ZHMIN ,  ZHMAX,  (DESC(l),  ..., 

DESC (8) ) ,  NW,  NMATIN ,  MATIN (1),  AMSS(l), 
AMOMX(l),  AMOMY(l),  AMOMZ(l),  AIE(l),  AKE(l), 

. . . ,  MATIN (NMATIN) ,  AMSS (NMATIN) ,  AMOMX (NMATIN, 
AMOMY (NMATIN) ,  AMOMZ (NMATIN) ,  AIE (NMATIN),  AKE 
(NMATIN) 

where 

CASE  is  the  EPIC  case  number 

NHELE  is  the  number  of  elements  to  transfer 
to  HULL 

NHELER  is  the  number  of  elements  per  record  on 
Tape  15 

STRAIN  is  the  value  of  the  SAIL  variable 

XHMIN,  ...,  ZHMAX  are  min  and  max  coordinates 
for  the  projectile  being  transferred 

DESC(l),  ...,  DESC (8)  is  the  problem  descrip¬ 
tion  from  the  dump  tape 

NW  is  the  number  of  words  per  element  being 
transferred  (27  if  STRAIN=0,  33  if  STRAIN=1) 

NMATIN  is  the  number  of  separate  materials 
being  transferred 

MATIN,  AMSS,  ...,  AKE  are  values  of  material 
number,  mass,  XYZ  linear  momenta,  total 
internal  energy  and  total  kinetic  energy  for 
all  of  the  materials  being  transferred 
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2, 


LAST 


Records,  each  of  which  contains  data  for 
NHELER  elements.  The  data  for  each  element 
is : 


MATERIAL  NUMBER 


U 

V 

W 


positions  for  the 
defining  nodes 


Element  centroid  velocity 
values  (i.e.,  one-fourth  the 
sum  of  the  values  at  the  four 
node  points) 


Mass  of  the  element 


Density  of  the  element 


Total  internal  energy  of  the  element 


Pressure 

X  and  Y  stress  deviators 
XY,  XZ,  YZ  shear  stresses 


ICHECK 

XYZ  normal  strains 

XY,  XZ,  YZ  shear  strains 


if  STRAIN  >  0 


e  -  accumulated  plastic  strain 
P 

X,  Y,  Z  values  for  all  EPIC 3  element  time 
history  points  in  the  region  being  transferred 
to  HULL 


User  input  to  HULLIN  is  seen  below: 
PROBLEM  DESCRIPTION 


CASE,  CYCLE 


MATERIAL  TRANSFER  DESCRIPTION 
NMATIN ,  MATIN (1) ,  . . . ,  MATIN (NMATIN) 


X,  Y,  Z  LIMITS  FOR  TRANSFER 

XMIN ,  XMAX,  YMIN ,  YMAX ,  ZMIN,  ZMAX 
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NMATIN 

is  the  number  of  materials  to  transfer 

MATIN (I) 

is  an  array  containing  material  numbers 
for  transfer 

XMIN ,  . . . ,  ZMAX 

define  spatial  regions  for  transfer 

To  facilitate  an  analysis  of  the  adequacy  of  the  transfer, 
program  HULLIN  prints  mass,  momenta  and  energies,  for  each  re¬ 
spective  material  being  transferred.  These  can  be  compared  to 
similar  data  provided  by  KEEL  as  it  processes  Tape  15. 
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SECTION  VIII 

CHANGES  TO  THE  STATE  PLOTTER 


The  state  plotter,  POST2 ,  was  minimally  changed  to  accomo¬ 
date  SAIL  options,  automatic  restart  tape  attachment  and  catalog¬ 
ing  of  Tape  99,  the  new  time  history  internal  laods  format  and  to 
provide  room  for  future  contour  variable  expansion.  In  addition, 
a  sizing  program,  P1SIZE,  was  written  to  read  SAIL  options  from 
restart  (or  time  history)  tapes,  check  LCM  availability  and  pro¬ 
vide  an  INPUT2  file  for  SAIL  processing  of  P0ST1. 

A  few  changes  were  made  in  input  in  addition  to  field  free 
format.  Case  number  must  be  read  as  the  first  input.  This  is 
required  in  order  to  attach  the  restart  tape.  Titles  for  plots 
have  been  dropped  and  the  problem  description  from  tape  inserted 
at  the  top  of  each  plot.  Internal  loads  are  now  specified  by  a 
card  beginning  with  a  100  (instead  of  15) .  This  change  was  made 
to  leave  room  for  future  expansion  of  the  contour  plots  variable 
array.  Of  course,  the  previously  terminating  blank  card  is  now 
a  card  with  a  single  zero  on  it. 

Input  for  POSTl  is  as  follows : 

CASE  NUMBER  CARD 


CASE 


GEOMETRY  PLOT  CARDS  AS  REQUIRED 


CFOMFTRY  PLOT  CARD  _ _ _ - _ _ 

| 

1,  VIEW,  CYCLE  IAXES ,  IFAIL 

PLOT  LIMIT  CARD  IF  IAXES=1  - - - - 

| 

XMAX,  XMIN ,  YMAX ,  YMIN ,  ZMAX,  ZMIN 

| 

3D  PERSPECTIVE  CARD  IF  VIEW=4  - 

XEYE ,  YEYE ,  ZEYE,  XPLANE ,  YPLANE ,  Z PLANE ,  LHIDE 

1 
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VELOCITY  VECTOR  PLOT  CARDS  AS  REQUIRED 


PLOT  CARD 


2,  VIEW,  CYCLE,  IAXES,  EDGE,  AROW, 

VS C ALE 

LIMIT 

CARD  IF  IAXES=1 

XMAX, 

XMIN,  YMAX,  YMIN,  ZMAX,  ZMIN 

3D  PERSPECTIVE  CARD  FOR  VIEW=4 

XEYE , 

YEYE ,  ZEYE,  XPLANE ,  YPLANE , 

Z PLANE 

CONTOUR  PLOT  CARDS  AS  REQUIRED _ 

PLOT  CARD _ 

TYPE,  VIEW,  CYCLE,  IAXES,  EDGE,  NLINE ,  PRINT,  ISYM 


LIMIT  CARD  IF  IAXES=1 _ 

XMAX ,  XMIN ,  YMAX ,  YMIN ,  ZMAX ,  ZMIN  ~1 

CONTOUR  SPECIFICATION  CARD _ 

VAR ( 1 ) ,  VAR ( 2 ) ,  ...,  VAR (NLINE) 


GEOMETRY,  VELOCITY,  CONTOUR  PLOTS  TERMINATING  CARD 


* 

INTERNAL  LOADS  PLOT  CARD(S)  AS  REQUIRED 


PLOT  CARD 


100, 

CYCLE,  IAXES,  IP, 

IV,  IM 

LIMIT 

CARD  FOR  IAXES=1 

PMAX , 

PMIN ,  VMAX ,  BMAX , 

BMIN ,  LENGTH 

INTERNAL  PLOTS  TERMINATING  CARD 


0 
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VIEW 

is  a  code  specifying  the  view  requested. 
VIEW=1,2,3  requests  two-dimensional 
plots  of  the  XZ,  YZ  and  XY  axes  respec¬ 
tively 

VIEW=4  requests  a  three-dimensional 
plot 

CYCLE 

is  the  cycle  number  of  the  dump  desired 
to  be  plotted 

IAXES 

determines  if  automatic  scaling  will 
prevail 

IAXES=0  means  automatic  scaling 

IAXES=1  means  that  coordinate  limits 
will  be  supplied 

IFAIL 

places  an  "X"  at  the  center  of  each 
failed  element  for  VIEW=1,2  or  3;  i.e., 
each  element  for  which  ICHECK  is  greater 
than  zero.  (ICHECK=-1  elements  are 

considered  as  separated  from  the  pro¬ 
blem  and  never  checked) 

XMAX ,  XMIN ,  YMAX , 

YMIN ,  ZMAX ,  ZMIN 

are  coordinate  limits  for  the  plot  and 
are  used  only  if  IAXES=1 

IEYE ,  YEYE ,  ZEYE 

are  coordinates  of  the  Observer 

XPLANE,  YPLANE ,  Z PLANE 

are  coordinates  included  in  the  plane 
for  which  results  should  be  plotted. 

The  plane  is  normal  to  a  line  from 
(XEYE ,  YEYE,  ZEYE)  to  (XPLANE,  YPLANE, 

Z PLANE) 

LHIDE 

specifies  if  lines  should  be  hidden 

LHIDE=0  means  plot  all  free  surfaces 

LHIDE=1  means  plot  only  those  surfaces  not 
hidden  from  the  observer 

(It  is  recommended  that  this  option  be  em¬ 
ployed  only  with  all  nodes  in  core  to 
avoid  excessive  CP  time) 
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EDGE 

specifies  whether  an  outline  should  be 
drawn  around  the  external  surfaces  for  a 
velocity  vector  plot 

EDGE=1  means  plot  the  outline 

EDGE=0  means  do  not  plot 

AROW 

if  set  to  1  will  draw  arrow  heads  on  the 
velocity  vectors 

VS C ALE 

is  a  multiplier  for  velocity  components. 

If  set  at  zero,  VSCALE  will  give  the 
longest  vector  a  length  of  two  percent  of 
the  vertical  axis 

TYPE 

specifies  the  type  of  contour  plot  desired 
TYPE  CONTOUR 

3  Pressure 

4-6  XYZ  normal  stresses 

7  XY  shear  stress 

8  XZ  shear  stress 

9  YZ  shear  stress 

10  Effective  stress  (VonMises  limit 

if  element  is  plastic) 

11  Equivalent  plastic  strain 

12  Temperature 

13  Internal  energy  per  unit  volume 

14  Plastic  work  per  unit  volume 

NLINE 

specifies  the  number  of  contours  to  be 
plotted.  This  number  is  currently  limited 
to  8.  If  NLINE=0,  a  default  option  is  in¬ 
voked  for  6  contours  at  values  of  5,  20, 

40,  60,  80  and  95  percent  of  the  range 
between  the  minimum  and  maximum  variable 
quantity  limits. 

PRINT 

prints  nodal  quantities  for  the  specified 
variable  if  set  to  1.  If  the  print  is  not 
desired  set  it  to  0. 

ISYM 

specifies  the  increment  at  which  symbols 
are  placed  on  contour  lines.  ISYM=1 
places  symbols  at  the  forward  end  of  each 
contour  line  within  an  element.  ISYM=5 
places  symbols  at  the  forward  end  of  each 
fifth  element,  etc. 

ISYM=0  will  not  put  any  symbols  on  the 
contour  lines. 
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VAR (I) 

specifies  the  magnitude  of  contours  to 
be  plotted  if  NLINE  is  greater  than  zero 

IP 

=  1  gives  plots  of  axial  loads 

IV 

=  1  gives  plots  of  shear  loads 

IM 

=  1  gives  plots  of  bending  moments 

PMAX,  PMIN 

are  bounding  coordinates  for  the  vertical 
axial  load  axis  (force) 

VMAX,  VMIN 

are  bounding  coordinates  for  the  vertical 
shear  load  axis  (force) 

BMAX,  BMIN 

are  bounding  coordinates  for  the  vertical 
bending  moment  axis  (force-distance) 

LENGTH 

is  the  maximum  coordinate  of  the  horizon¬ 
tal  axis  which  is  defined  as  the  deformed 
centerline  distance  from  the  free  end  of 
the  rod  (rod  length) 
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SECTION  IX 

TIME  HISTORY  PLOTTING 


The  time  history  plotter,  POST2 ,  as  well  as  the  method  of 
saving  time  history  data  were  all  but  completely  rewritten.  Some 
P0ST2  rewrite  was  required  because  of  the  method  of  packing  data 
on  the  time  history  file.  However,  a  large  amount  of  rewrite  was 
required  because  of  the  old  P0ST2  requirement  that  data  for  four 
types  of  time  history  stations  had  to  fit  in  addressable  core. 

This  imposed  severe  restrictions  on  the  amount  of  data  which  could 
be  plotted  during  any  single  run.  Titles  in  A6  format  also  had 
to  be  changed  to  provide  reasonable  labels  on  a  60  bit  word  size 
machine. 

The  approach  taken  to  processing  the  data  was  to  employ 
separate  files  for  storage  of  element,  node,  system  and  internal 
loads  data.  In  particular,  when  data  for  an  element,  a  node, 
the  system  or  a  plane  for  internal  loads  is  requested,  POST2  reads 
the  EPIC  time  history  tape  and  places  all  data  for  that  element, 
node,  etc.  on  a  file.  This  file  is  then  read  by  the  appropriate 
plotting  subroutine  and  the  desired  variable  (pressure,  stress 
components,  etc.)  is  plotted.  Because  of  this  processing  techni¬ 
que,  it  is  desirable  to  call  for  all  required  plots  for  a  parti¬ 
cular  element  node  or  plane  prior  to  switching  the  request  to  a 
new  element,  node  or  plane. 

The  requirement  to  submit  a  plot  title  with  each  plot  request 
has  been  eliminated.  Instead,  the  problem  description  assigned  in 
PREP  will  automatically  be  recovered  from  the  time  history  tape 
and  placed  at  the  top  of  each  plot. 

The  addition  of  strains  to  data  carried  in  elements  (if  the 
SAIL  option  STRAIN  is  greater  than  zero)  allowed  us  to  change 
P0ST2  to  offer  time  history  plots  of  exx,  eyy,  ezz,  exy,  exz,  eyz 

and  the  maximum  principle  strain  for  elements.  In  addition,  if 
PYFAIL  is  greater  than  zero  in  the  problem  calculation,  P/Y  can 
be  plotted  as  a  function  of  time  for  an  element. 

The  P2SIZE  program  is  run  prior  to  running  POST2 .  This  pro¬ 
gram  reads  the  SAIL  option  block  from  the  time  history  tape  and 
prepares  an  INPUT2  file  for  SAIL  to  use  when  processing  POST2 . 

Similarly  to  PREP,  MAIN  and  POST1,  POST2  will  automatically 
attach  the  latest  time  history  file  if  CATPLT  is  greater  than 
zero.  In  addition,  it  will  automatically  catalog  the  Tape  99 
file  if  CATPLOT  is  greater  than  zero.  The  file  is  cataloged  on  the 
A  machine  at  BRL  under  a  directory  specified  by  the  value  of  the 
SAIL  option  USER.  The  file  is  cataloged  as: 

P  2  XXXXPTYYYY 
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where  the  problem  (case)  number  is  XXXX.YYYY.  If  the  time 
history  tape  is  a  real  tape  it  is  recommended  that  it  be  copied 
to  a  disc  in  the  control  stream  because  of  the  many  potential 
READ'S  and  REWIND'S  which  can  occur  in  P0ST2. 


Input  for  the  new  P0ST2  is  as  follows: 
CASE  NUMBER  CARD _  _ 


CASE 

SYSTEM  PLOT 

CARD ( S )  AS  REQUIRED 

TYPE,  IAXES , 

SCALE,  TMAX , 

TMINi 

VMAX , 

VMIN 

INTERNAL  LOAD  CARD(S)  AS 

REQUIRED 

TYPE,  IAXES, 

LAYER,  SCALE 

,  TMAX 

,  TMIN 

,  VMAX 

,  VMIN 

NODE  PLOT  CARDS  AS  REQUIRED 

TYPE,  IAXES, 

NODE ,  SCALE , 

TMAX, 

TMIN , 

VMAX, 

VMIN 

ELEMENT  PLOT 

CARDS  AS  REQUIRED 

TYPE,  IAXES, 

ELE,  SCALE, 

TMAX,  TMIN , 

VMAX, 

VMIN 

ENDING  CARD 

0 

TMAX,  TMIN,  VMAX,  VMIN 

these  min  and  max  values  for  each  plot  will 
be  used  if  IAXES=1.  If  IAXES=0,  min  and 
max  values  are  automatically  chosen.  Note 
that  values  for  TMAX,  TMIN,  VMAX,  and  VMIN 
must  be  present  on  the  card  even  if  IAXES=0. 
In  this  case,  any  values  -  including  zeros 
are  acceptable. 

SCALE 

is  a  multiplier  on  the  variable  specified 
by  the  input  card 

NODE 

is  specific  node  number 

ELE 

is  specific  element  number 

LAYER 

is  specific  layer  or  plane  number  for 
internal  loads  plots 

TYPE 

specifies  the  variable  to  be  plotted 
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The  acceptable  values  are: 
TYPE 


MEANING 


1 

SYSTEM 

Total  Energy 

2 

SYSTEM 

Kinetic  Energy 

3 

A 

SYSTEM 

Internal  Energy 

4 

SYSTEM 

Plastic  Energy 

5-7 

SYSTEM 

XYZ  Maximum  Coordinates 

8-10 

SYSTEM 

XYZ  Minimum  Coordinates 

11-13 

SYSTEM 

XYZ  Centers  of  Gravity 

14-16 

SYSTEM 

XYZ  Linear  Momenta 

17-19 

SYSTEM 

XYZ  Velocities 

20 

SYSTEM 

Total  Velocity 

21-23 

SYSTEM 

XYZ  Angular  Momenta 

24-26 

SYSTEM 

XYZ  Angular  Velocities 

27 

SYSTEM 

Axial  Load 

28 

SYSTEM 

Shear  Load 

29 

SYSTEM 

Bending  Load 

30-32 

NODE 

XYZ  Positions 

33-35 

NODE 

XYZ  Velocities 

36-38 

NODE 

XYZ  Accelerations 

39 

NODE 

Node  Pressure 

40 

ELEMENT 

Pressure 

41-43 

ELEMENT 

XYZ  Normal  Stresses 

44 

ELEMENT 

XY  Shear  Stress 

45 

ELEMENT 

XZ  Shear  Stress 

46 

ELEMENT 

YZ  Shear  Stress 

47 

ELEMENT 

Effective  Stress  (Von  Mises 
limit  if  element  in  plastic 
flow) 

48 

ELEMENT 

ep-equivalent  plastic  strain 

49 

ELEMENT 

Temperature 

50 

ELEMENT 

Internal  Energy  per  unit 
volume 

51 

ELEMENT 

Plastic  work  per  unit 
volume 

52-54 

ELEMENT 

XYZ  Normal  Strains 

55 

ELEMENT 

XY  Shear  Strain 

56 

ELEMENT 

XZ  Shear  Strain 

57 

ELEMENT 

YZ  Shear  Strain 

5  8 

ELEMENT 

Maximum  Principle  Strain 

59 

ELEMENT 

P/Y 

NOTE: 

These  last  two  quantities  are 
run  was  made  with  IPNODE=l. 

nodal  averages  if  the  MAIN 
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SECTION  X 

EPIC  MONITOR  SYSTEM 


Routines  were  developed  to  allow  the  user  to  communicate  with 
a  running  EPIC3  job  for  job  status  determination  and  job  control. 
The  routines  were  patterned  after  those  in  use  in  the  HULL  system. 
Communications  is  provided  through  a  file  catalogued  as  EPIC3C0N- 
TROL  on  sysset  on  the  BRL  Z  machine.  The  file  is  catalogued  under 
the  directory  KIMSEY. 

Communications  is  currently  limited  to  program  MAIN,  but  it 
would  be  a  simple  matter  to  add  the  necessary  routines  to  PREP  or 
the  plotting  programs.  Subroutine  CONTW  in  MAIN  (actually  brought 
into  MAIN  from  PROLOGUE)  is  called  at  the  end  of  each  cycle  if  the 
SAIL  option  RCONTRL  is  greater  than  zero.  Subroutine  CONTW  checks 
that  the  file  EPIC3C0NTR0L  exists  by  trying  to  attach  it.  If  the 
file  does  not  exist,  it  creates  one.  The  random  access  file  is 
created  using  mass  store  commands  -  so  that  attaches  and  catalogs 
can  be  made  without  accompanying  system  printouts. 

The  file  is  currently  configured  for  twenty-one  fifty  word 
records.  The  first  record  is  a  table  containing  case  numbers  for 
the  following  data  records.  For  example,  if  the  fifth  word  in  the 
first  record  is  2.5,  then  the  fifth  data  record  (sixth  record  on 
the  tape)  contains  the  data  for  problem  (case)  2.5.  A  zero  in 
the  first  twenty  words  of  the  first  record  indicates  an  available 
data  record. 

Each  data  record  for  a  running  (or  previously  run)  problem 
contains  the  following  information: 

WORD  NUMBER  DATA 


1 

Current  value  of  NCYCLE 

2 

Current  value  of  TIME 

3 

Current  value  of  DT 

4 

Controlling  element,  LCRIT 

5-12 

Problem  description 

13 

Central  Memory  being  used 

14 

LCM  being  used 

15 

CP  time  used  during  the  current  run 

16 

CP  time  per  element  per  cycle  for  the  run 
this  time 

to 

17 

Date  of  last  file  update  for  the  run 

18 

Time  of  day  for  last  file  update 

19 

The  number  of  elements  in  the  problem,  NELE 

20 

The  number  of  nodes  in  the  problem,  NNODE 

21 

The  number  of  system  data  records  written 
the  time  history  tape 

on 

22 

The  number  of  internal  loads  records 
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The  number  of  node  records 
The  number  of  element  records 
The  value  of  LCMNODE 
The  value  of  LCMNBUF 
The  value  of  LCMELEM 
The  value  of  LCMEBUF 

The  number  of  node  blocks  in  the  problem,  NNBLK 
The  number  of  element  blocks  in  the  problem,  NEBLK 
The  current  value  of  NBAND 
The  current  value  of  LBAND 
Not  used 
NHULL 

MODE=EPIC  or  END  depending  on  whether  the  problem 
has  terminated  or  not 
Not  used 

Control  word=0  if  no  control  has  yet  been  exercised 
=1  if  the  problem  should  be  terminated 
=2  if  a  printout  of  the  output  file  to  this  point 
is  desired 

(the  control  word  is  input  from  program  EMON  which 
the  user  must  interactively  run  to  monitor  EPIC3 
runs) 

As  can  be  seen,  there  is  sufficient  room  for  a  great  deal  more 
monitoring  information  and/or  control  options.  Subroutine  CONTW 
and  program  EMON  are  sufficiently  commented  so  that  they  could 
easily  be  modified  by  the  BRL  to  provide  more  data  or  control  as 
desired. 

Program  EMON  is  a  separate  program  on  the  EPIC  SYSTEM  file. 

It  can  be  separated  from  the  file  by  running  the  following  control 
cards  on  the  A  machine. 

Job  card 

Account  card 

ATTACH,  HULLIB,  ID=KIMSEY. 

LIBRARY,  HULLIB. 

ATTACH,  OLD,  EPICSYS,  ID=KIMSEY. 

COPYCR,  INPUT,  A. 

REWIND,  A. 

DYTHUL,  I=A. 

CATALOG,  SAIL,  EMON,  ID=XXX. 

*EOR 

SAIL  PROGRAM  EMON 

As  a  result  of  the  run,  the  source  file  for  EMON  will  be  catalogued 
as  file  EMON  on  sysset  directory  XXX.  The  file  can  then  be  at¬ 
tached  from  a  terminal  and  compiled  on  the  A  machine  using  the 
STATIC  option.  It  must  be  run  with  the  HULLIB  library  attached 
and  declared.  The  Interactive  compile  command  to  invoke  STATIC  is: 

FTN ,  STATIC,  I=EMON,  B=RUN 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36-49 

50 
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where  RUN  is  the  user's  desired  name  for  the  load  and  execute  file 
(the  default  is  LGO) . 

Program  EMON  connects  INPUT  (Tape  5)  and  OUTPUT  (Tape  6)  to 
the  machine  on  which  it  is  running.  Thus  running  interactively 
the  user  inputs  requests  on  a  remote  terminal  and  the  program 
output  is  automatically  printed  on  the  terminal. 

Commands  and  data  displayed  from  EMON  are  as  illustrated 
below: 

COMMAND: LIST 

EMON:  EMON  prints  a  list  of  EPIC  problems  on  the  EPIC3C0NTR0L 

file.  For  each  problem,  EMON  prints  status  (MODE) , 
cycle,  problem  time  and  time  increment. 

COMMAND: DISPLAY  XXXX.YYYY 

EMON:  EMON  prints  the  data  in  the  problem  data  record  for 

case  XXXX.YYYY  displaying  problem  description,  cycle, 
time,  time  increment,  NHULL,  NELE,  NNODE,  etc. 

COMMAND: DIRECT  XXXX.YYYY  PARTOUT 

EMON:  EMON  places  a  2  in  the  control  word  for  case  XXXX.YYYY 

and  when  EPIC3  calls  CONTW  it  will  use  this  control 
word  to  return  the  output  file  making  it  available  to 
the  user. 

COMMAND: DIRECT  XXXX.YYYY  STOP 

EMON:  EMON  places  a  1  in  the  control  word  for  case  XXXX.YYYY 

and  when  EPIC3  calls  CONTW  it  uses  this  control  word 
to  cause  EPIC  termination. 

COMMAND : STOP 

EMON:  EMON  ceases  execution. 
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SECTION  XI 

SAMPLE  LINKED  PROBLEM 


✓ 


A  sample  linked  HULL/EPIC  problem  was  run  to  demonstrate  pro¬ 
per  operation  of  the  linking  algorithms.  The  calculation  with  a 
comparison  to  experimental  results  will  be  presented  in  a  separate 
technical  report.  In  this  section  we  will  primarily  discuss  the 
input  decks  used  for  the  EPIC  portion  of  the  problem. 

Figure  1  presents  density  contours  from  the  HULL  calculation 
of  an  L/D=10  penetrator  impacting  a  thin  armor  plate  at  approxi¬ 
mately  1  km/sec  velocity.  Figure  2  presents  EPIC  velocity  vec¬ 
tors  at  various  times  in  the  EPIC  portion  of  the  calculation. 
Forward  element  planes  in  the  EPIC  calculation  were  dropped  at 
36.5  and  42  microseconds  resulting  in  the  shortened  rod  seen  at 
60  microseconds.  The  planes  were  dropped  because  they  became  too 
highly  compressed  and  were  essentially  stopping  the  calculation. 

The  initial  HULL  and  EPIC  geometries  are  seen  in  Figure  3 . 
Since  the  hemispherical  nose  cap  of  the  rod  is  eroded  away  early 
in  the  HULL  calculation  it  was  not  modelled  at  all  in  the  EPIC 
calculation.  This  does  not  affect  the  accuracy  of  the  calcula¬ 
tion  since  HULL  velocities  used  to  drive  EPIC  will  be  inside  the 
rod  at  the  correct  locations  and  since  the  cap  area  will  have  been 
removed  from  the  penetrator  prior  to  termination  of  the  HULL  cal¬ 
culation.  There  were  541  HULL  stations  in  the  rod  nose  area. 

The  rod  modelled  in  EPIC  has  840  nodes  and  3132  elements. 

The  EPIC  calculation  required  2370  CP  seconds  to  reach  a  problem 
time  of  200  microseconds  in  1530  cycles.  Overall,  the  EPIC 
running  time  on  the  CDC7600  at  BRL  can  be  computed  to  be  1.84 
billiseconds  per  node  per  cycle  (all  of  the  nodes  were  contained 
in  LCM  for  the  calculation) .  This  is  slightly  greater  than  the 
usual  1.4  milliseconds  because  of  the  extra  station  location 
calculations  required  because  of  the  link. 

Figures  4  and  5  present  PREP  and  MAIN  input  decks  for  the 
calculation  in  CGS  units.  As  seen  in  these  figures  the  rod  nose 
was  located  at  0.385  cm  (since  the  nose  cap  was  not  included)  and 
the  rod  terminated  at  a  Z  value  of  7.7  cm.  It  has  a  constant 
radius  of  0.385  cm.  The  staballoy  was  modelled  using  the  HULL 
equation  of  state  with  a  P/Y  failure  curve  determined  from  notched 
tensile  tests.  The  SAIL  variable  FAIL  was  set  to  zero  for  the  cal¬ 
culation.  Most  Of  the  MAIN  deck  consists  of  EPIC  time  history 
station  location  data  (34  element  stations  in  all) .  The  deck 
specifies  that  a  dump  be  made  every  20  microseconds  and  that  the 
run  be  terminated  at  200  microseconds  or  at  the  time  when  the  min¬ 
imum  Z  value  of  the  rod  minus  1.E5  X  TIME  reaches  -20.085  cm. 

This  is  the  location  of  the  second  target  plate.  Elements  numbered 
3025  and  greater  (the  forward  rod  plane)  will  be  dropped  at  36.5 
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T  =  41  ys 


T  =  47  ys 


igure  1 . 


OTI/BRL  HULL  Calculation  of  RHA  Perforated  by 
staballoy  at  a  velocity  of  1  km/sec  and  an  lm 
pact  angle  of  65  degrees. 
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Figure  2.  EPIC3  velocity  Plots  for  Linked  Calculation. 
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Figure  4 .  PREP  Input  Data  for  the  Sample  Problem 
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Figure  5.  MAIN  Input  Data  for  the  Sample  Problem 
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microseconds.  Elements  numbered  2917  and  greater  (the  second 
most  forward  plane)  will  be  dropped  at  42  microseconds.  The 
automatic  element  drop  routine  which  reacts  to  element  time  step 
was  not  used. 
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